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THE GLOBAL NONLINEAR STABILITY OF MINKOWSKI SPACE 
FOR SELF-GRAVITATING MASSIVE FIELDS. 

The Wave-Klein-Gordon Model 

PHILIPPE G. LEFLOCH AND YUE MA 


Abstract. The Hyperboloidal Foliation Method (introduced by the authors in 2014) is ex¬ 
tended here and applied to the Einstein equations of general relativity. Specifically, we establish 
the nonlinear stability of Minkowski spacetime for self-gravitating massive scalar fields, while 
existing methods only apply to massless scalar fields. First of all, by analyzing the structure of 
the Einstein equations in wave coordinates, we exhibit a nonlinear wave-Klein-Gordon model 
defined on a curved background, which is the focus of the present paper. For this model, we 
prove here the existence of global-in-time solutions to the Cauchy problem, when the initial data 
have sufficiently small Sobolev norms. A major difficulty comes from the fact that the class of 
conformal Killing fields of Minkowski space is significantly reduced in presence of a massive 
scalar field, since the scaling vector field is not conformal Killing for the Klein-Gordon operator. 
Our method relies on the foliation (of the interior of the light cone) of Minkowski spacetime by 
hyperboloidal hypersurfaces and uses Lorentz-invariant energy norms. We introduce a frame 
of vector fields adapted to the hyperboloidal foliation and we establish several key properties: 
Sobolev and Hardy-type inequalities on hyperboloids, as well as sup-norm estimates which cor¬ 
respond to the sharp time decay for the wave and the Klein-Gordon equations. These estimates 
allow us to control interaction terms associated with the curved geometry and the massive field, 
by distinguishing between two levels of regularity and energy growth and by a successive use of 
our key estimates in order to close a bootstrap argument. 
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1. Introduction 

1.1. The global existence problem. In this paper and its companion [53], we study the global- 
in-time existence problem for small amplitude solutions to nonlinear wave equations, with a two¬ 
fold objective: 

• First, we provide a significant extension of the Hyperboloidal Foliation Method, recently 
proposed by the authors [55]. This method is based on a foliation of the interior of the 
future light cone by hyperboloidal hypersurfaces and on Sobolev and Hardy inequalities 
adapted to this foliation. This method takes its root in work by Klainerman m and, 
later on, Hdrmander m concerning the standard Klein-Gordon equation. In comparison 


^Published in: Communications in Mathematical Physics (2016). 
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to our earlier theory in [22] , we are now able to encompass a much broader class of coupled 
wave-Klein-Gordon systems. 

• Our second objective is to apply this method to the Einstein equations of general relativity 
and arrive at a new approach for proving the nonlinear stability of Minkowski spacetime. 
Our method covers self-gravitating massive scalar fields (as will be presented in full details 
in [28)1. while earlier works were restricted to vacuum spacetimes or to spacetimes with 
massless scalar fields; cf. Christodoulou and Klainerman [7], and Lindblad and Rodnianski 
PSI [?7] , as well as Bieri and Zipser [3] . 

The problem of the global dynamics of self-gravitating massive fields had remained open until now. 
The presence of a mass term poses a major challenge in order to establish a global existence theory 
for the Einstein equations (and construct future geodesically complete spacetimes). Namely, the 
class of conformal Killing fields of Minkowski spacetime is reduced in presence of a massive scalar 
field, since the so-called scaling vector field is no longer conformal Killing and, therefore, cannot 
be used in implementing Klainerman’s vector field method mm- 

In suitably chosen coordinates, the Einstein equations take the form of a coupled system of 
nonlinear wave-Klein-Gordon equations. More precisely, as in [27] . we introduce wave coordinates, 
also called harmonic or De Bonder gauge [1], which allows one to exhibit the (quasi-null, see below) 
structure of the Einstein equations. The Hyperboloidal Foliation Method [H] was introduced 
precisely to handle such systems. Yet, due to the presence of metric-related terms in the system 
under consideration, an important generalization is required before we can tackle the Einstein 
equations. Proposing such a generalization is our main purpose in the present paper. 

By imposing asymptotically flat initial data on a spacelike hypersurface with sufficiently small 
ADM mass, one can first solve the Cauchy problem for the Einstein equations within a neigh¬ 
borhood of this hypersurface (see [12] for a sketch of the argumenlQ) and, next, formulate the 
Cauchy problem when the initial data are posed on a hyperboloidal hyperspace or, alternatively, 
on a hyperboloid for the flat Minkowski metric after introducing suitable coordinates. In fact, 
the hyperboloidal Cauchy problem is, both, geometrically and physically natural. More precisely, 
let us consider Minkowski spacetime in standard Cartesian coordinates and observe 

that points on a hypersurface of constant time t cannot be connected by a timelike curve, while 
points on a hyperboloid can be connected by such curves. Hence, hyperboloidal initial data can 
be “physically prepared”, while data on standard flat hypersurfaces cannot. An alternative stand¬ 
point would be to pose the Cauchy problem on a light cone, but while it is physically appealing 
and the Cauchy problem on a light cone has not been proven to be convenient for global analysis. 

We emphasize that hyperboloidal foliations were used by Friedrich HU da in order to establish 
the stability of Minkowski space in the future of a hyperboloidal hypersurface. Hyperboloidal 
foliations have also been found to be very efficient in numerical computations [TOllllllSOllM]. 

As was demonstrated in [22] for a rather general class of nonlinear wave equations, analyzing 
the global existence problem is quite natural in the hyperboloidal foliation of Minkowski spacetime 
and, importantly, lead to uniform bounds on the energy of the solutions. Before proceeding with 
further details, let us summarize the main features of the method we propose: 

• Lorentz invariance. We rely on the foliation of Minkowski space by hyperboloids (de¬ 
fined as the level sets of constant Lorentzian distance from some origin), so that the 
fundamental energy of the wave-Klein-Gordon equations remains invariant under Lorentz 
transformations of Minkowski spacetime. Observe that in our construction, all the hyper¬ 
boloids are asymptotic to the same limiting cone and approach the same sphere at infinity. 
(In particular, no energy can escape through null infinity.) 

• Smaller set of Killing fields. We avoid using the scaling vector field S := rdr + tdt, 
which is the key in order to handle Klein-Gordon equations and cover the Einstein-matter 
system when the evolution equation for the matter is not conformally invariant. 


^The time of existence of the solution can be made arbitrarily large for compactly supported initial data with 
sufficiently small norm, so that this neighborhood does contain a hyperboloidal hypersurface. 
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• Sharp rate of time decay. In order to control source-terms related to the curved 
geometry, we establish sharp pointwise bounds for solutions to wave equations and Klein- 
Gordon equations with source-terms. 

In the rest of this introduction, we explain how to derive, from the Einstein equations, a model 
problem which will be our main focus in the present paper. 

1.2. Einstein equations for massive scalar fields. We thus consider the Einstein equations 
for an unknown spacetime (M, g): 

R 

(1-1) f^a/3 TafUj 

where Ra /3 denote^ the Ricci curvature tensor and R = g^^^Rap denotes the scalar curvature. 
The matter is taken to be a massive scalar field with potential V = V {(j>) and stress-energy tensor 

( 1 . 2 ) := 

and, specifically, 

(1.3) V{^) := 

where > 0 represents the mass of the scalar field. By applying V“ to (US and using the Bianchi 
identity 

V“(R„;3 - (R/2)g„;3) = 0, 

we easily check that the Einstein-scalar field system implies 
(1.4a) Ra/3 = Sir{ya4''^i34> + V{(j)) gap), 

(1.4b) = W(</)) = cV- 

The Cauchy problem for the Einstein-scalar field equations is formulated as follows [5]. An 
initial data set consists of a Riemannian three-manifold (M,g), a symmetric two-tensor K defined 
on M, and two scalar fields {(j>Q,(j>i) defined on M. We then seek for a (3 -t l)-dimensional 
Lorentzian manifold (M, g) satisfying the following properties: 

• There exists an embedding i : M ^ M such that the induced metric i*{g) coincides with g, 

while the second fundamental form of i{M) c M coincides with the prescribed two-tensor 
K. _ 

• The restriction of (j) and to i{M) coincides with the data 4>o and (j)i respectively, where 

denotes the (future-oriented) unit normal to i{M) c M. 

• Moreover, the manifold {M,g) satisfies the Einstein equations (11.41) . 

More precisely, one seeks for a globally hyperbolic development of the given initial data, that 
is, a Lorentzian manifold such that every time-like geodesic extends toward the past direction in 
order to meet the initial hypersurface M. Furthermore, a notion of maximal development was 
defined by Choquet-Bruhat and Geroch [61E] and such a development was shown to exist for a 
large class of matter models. The maximal development need not be future geodesically complete, 
and a main challenge in the field of mathematical general relativity is the construction of classes 
of future geodesically complete spacetimes. 

Furthermore, it should be emphasized that, in order to fulfill the equations (d, the initial 
data set (M,g, iL) cannot be arbitrary and must satisfy Einstein’s constraint equations: 

R + K,, - (Kl)^ = SttToo, 

( 1 - 5 ) , 

V = SttToj, 

where R is the scalar curvature of the metric g and V denotes its Levi-Givita connection, and the 
terms Tqo and T^i are determined from the data (j)^ and (j)^. 

^Throughout, Greek indices a, /3, 7 take values 0,1, 2, 3 and Einstein convention is used. 
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Minkowski spacetime provides one with a trivial solution to the Einstein equations, which 
satisfies the Cauchy problem associated with the initial data when M = is 

endowed with the standard Euclidian metric and K = 0, while the matter terms vanish identically 
4‘o ~ = 0. The question we address in the present paper is whether this solution is dynamically 

stable under small perturbations of the initial data. More precisely, given an initial data set 
(M, 5 , K, such that M is diffeomorphic to 'g is close to the flat metric and K, (f)^ are 

sufficiently small, does the associated solution {M,g) to the Einstein-massive scalar field system 
remain close to the flat Minkowski spacetime ? 

Clearly, this nonlinear stability problem is of fundamental importance in physics. It is expected 
that Minkowski spacetime is the ground state state of the theory with the lowest possible energy. 
As far as massless scalar fields are concerned, the nonlinear stability of Minkowski spacetime was 
indeed established in Christodoulou and Klainerman’s pioneering work [7]. In the present work 
(including [23)1. we solve this question for massive scalar fields. 


1.3. Einstein-scalar field equations in wave coordinates. Our first task is to express the 
field equations (11.41) in a well-chosen coordinate system and then derive our wave-Klein-Gordon 
model problem. We follow [3 [26] and work in wave coordinates satisfying, by definition, 

(1.6) □ga;“ = 0. 

We postpone to [23] the details of the derivation and directly write the formulation of the Einstein- 
massive scalar field equations in wave coordinates: 

(1.7a) ngga /3 = Qap{g\dg,dg) + Pap{g;dg,dg) - 167r(d„(/>d/3^ + E(0)5a^), 


(1.7b) (</>) = «, 

where Dg := g°‘^dadp is referred to as the (reduced) wave operator in curved space. In ()1.7)1 . we 
distinguish between several types of nonlinearity: 


( 1 . 8 ) 


(1.9) 


• Null terms. The quadratic terms 

Qap '■= g^^ g^^ Ssga\'ds'gp\ — g^^ g^^ {dsgaX’Sxgps’ — Ssgps’S\ga\') 

+ {dagx'S'dsgxp — dagxpdsgx'S') + -^g^^ g^^ {SagxpSygss' - dagss'dx'gxp) 

+ g^^ g^^ {Spgx'S'dsgxa — S/sgxa^sgx'S') + 2^^^ {^/sgxa^x'gss' — dpgss'dygxa) 

are standard null forms with cubic corrections. Their treatment in a global existence proof 
is a now classical matter and, in particular, are already dealt with by standard methods. 

• Quasi-null terms. The quadratic terms 

Pap ■■= g^^'dagsx'Spgxs' + ^g^^'g^^'dpgss'dagxx' 


are referred to as “weak null” terms in [26] , but we prefer to propose the new terminology 
“quasi-null terms”. As first noted in [26] . quasi-null terms are found to be analogous to 
standard null terms, provided the tensorial structure of the Einstein equations and the 
wave coordinate condition are carefully taken into account. 

• Curved metric terms. Setting now 


( 1 . 10 ) 




hap :— rUap gap 


and considering the term \Z\ggap, we see that we must also treat the quasi-linear terms 

^ da' dp> hap , ^ da' dp'cj). 


We will deal with these metric-related terms by the following two approaches: 

— Eirst, thanks to the wave coordinate condition, we can assume that behaves 
essentially like a null quadratic form and consider, therefore, that h°‘^ is null. More 
precisely, the first term 

^ da'dp'hap, h°‘^ being a null form 
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can be treated by our arguments in [22) . 

— The second quasi-linear term h°‘ ^ da'dpKp (without necessarily imposing the null 
condition) requires our new technique which is presented in this paper and is based on 
sharp sup-norm bounds for solutions to wave equations and Klein-Gordon equations. 

Our aim is presenting first in a simplifed form several arguments that will be required to cope 
with the full system of Einstein equations in [28] . In order to derive here a model problem, we 
proceed by removing (from (11.71) 1: 

• the null terms Qap (which are handled in [22)1. 

• the quasi-null terms Pa /3 (postponed to [23], where the structure of the Einstein equations 
and the wave coordinate condition will be discussed), and 

• the quasi-linear terms ^ da'dpihap (to be treated by the wave coordinate condition and, 
in turn, the method already presented in [22)1. 

These formal simplifications, therefore, lead us to the mode0 

nha/B = da(j)dp4> + rnal3V{(j)), 
n^ = H^^{h)dc,dpcj}+V'{cl}), 

with unknowns hap,4’ defined on Minkowski space, where can be assumed to depend 

linearly on hap- We are primarily interested in the potential and, therefore after 

changing the notation, we arrive at the following system of two coupled equations: 

— = P°‘^daVdi3V + 

— Du -f H°‘^udadpv + c^v = 0, 

where u, v are two scalar unknowns and P“^, R, c are given constants (and only the obvious 
positivity condition > 0 is relevant). 

1.4. Analysis on the model problem. As illustrated by the derivation above, in order to deal 
with the Einstein-massive scalar field equations, we must weaken a key assumption made in [22) 
and, as we will see, cope with wave equations posed on a curved space for which the Minkowski 
metric need not represent the underlying geometry in a sufficiently accurate manner. Namely, we 
recall that, in the notation of [221 Section 1], interaction terms like ui^ddvj involving components 
of wave equations and component vj of Klein-Gordon equations were not included in our theory. 
The same restriction was also assumed in a pioneering work by Katayama [Min] on wave-Klein- 
Gordon equations. In the present paper, we overcome this challenging difficulty and extend our 
earlier analysis (of the system (1.2.1) in [22] by now removing the condition (1.2.4e) therein). 

To this end, in the present paper, we derive and take advantage of two pointwise estimates: 

• A sharp sup-norm estimate for solutions to the wave equation in Minkowski 
space with source-term, as stated in Theorem 13.11 below. Suitable decay is assumed on 
the source-term, as is relevant for our analysis, and the proof is based on the solution 
formula available for the wave equation in flat space. 

• A sharp sup-norm estimate for solutions to the Klein-Gordon equation in curved 
space in (3-l-l)-dimensions (as stated in Theorem l3.3l below). Our estimate is motivated by 
a pioneering work by Klainerman [19] on the global existence problem for small amplitude 
solutions to nonlinear Klein-Gordon equations in four spacetime dimensions. 

Note that an estimate as above could also be derived in (24- l)-dimensions with different rates |5Sj. 
Global existence results for nonlinear Klein-Gordon equations were also established by Shatah in 
the pioneering work m- Klein-Gordon systems have received a lot of attention in the literature 
and we can, for instance, refer to [iisjiHi m in HI and the references therein. 


^Our convention for the wave operator is □ := —dtdt + 2a 
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For clarity in the presentation, we do not treat the most general class of systems, but based 
on our formal derivation from the Einstein equations, we now study the following wave-Klein- 
Gordon model: 

— Du = P°'^davdpv + Rv^, 

^ ^ + u H'^^dadpv + c^v = 0, 

with unknowns u, v posed on Minkowski space and prescribed initial dated uq, ui, vg, vi posed 
on the spacelike hypersurface t = 2: 


( 1 . 12 ) 


U\t=2 = Uq, dtu\t=2 = Ui, 

V\t =2 = Vo, dtv\t =2 = Vl. 


Here, P“^, R, c are given constants, and the initial data are sufficiently smooth functions that 
are compactly supported in the unit ball {(xi)^ + {x 2 )‘^ + (xs)^ < 1} with x = {xi,X 2 ,X 3 ) e 
We emphasize that, according to our analysis in Section ITT^ (ll.ll|) includes the essential dif¬ 
ficulty arising in the Einstein-massive field system. Note in passing that, in (11.111) . there is no 
such term like Ru^ which would imply finite time blow-up (as first pointed out by John [15)1. Our 
main result in the present paper is as follows. 


Theorem 1.1 (Global existence theory for the wave-Klein-Gordon model). Consider the nonlinear 
wave-Klein-Gordon system dnii for some given parameter values P“^, R, H°‘^ and c > 0. Given 
any integer N ^ 8, there exists a positive constant eg = £o{N) > 0 such that if the initial data 
satisfy 

(1-13) ||(Mo,'yo)||H«‘+i(B3) + II(mi,'Ci)||hn-(r3) < £0, 

then the Cauchy problem (I1.11|) - (I1.12I1 admits a global-in-time solution. 

As done in [23) . the Gauchy problem can be reformulated with initial data prescribed on a 
hyperboloid and the smallness condition (I1.13|) leads to a similar smallness condition for the 
hyperboloidal initial data. As already pointed out in [22], the presence of the quasi-linear term 
uH°‘^do,dpv may possibly change the asymptotic behavior of solutions for large times. In fact, our 
proof will only show that the lower-order energy of the wave component remains globally bounded 
for all times, while the high-order energy of the wave component u and the lower-order energy 
of the Klein-Gordon component v could in principle grow at the rate t^ for some (small) d > 0. 
On the other hand, the higher-order energy associated with the Klein-Gordon component v may 
significantly increase at the rate t^+P'^ for some (small) d > 0. 

The proof of Theorem 1 1.1 1 will occupy the rest of this paper which we outline as follows: 

• Proceeding with a bootstrap argument, we assume that, within some hyperbolic time 
interval, the hyperboloidal energy of suitable derivatives of the unknowns (up to a certain 
order) satisfy a set of bounds. 

• Our assumptions use two levels of regularity and distinguish between the behavior of 
lower-order and higher-order energy norms, the low-order derivatives enjoying a much 
better control in time. Recall that, in [35], we could already prove that the lower-order 
energy of the wave component is uniformly bounded in time, but the growth rate for the 
high-order Klein Gordon energy was solely t^. 

• By Sobolev inequality (on hyperboloids), we can turn these type inequalities to a set 
of sup-norm estimates, which we refer to as basic decay estimates. These decay estimates 
are not sharp enough in order to close our bootstrap argument. 

• Relying on these basic decay estimates, we establish refined decay estimates by relying on 
two technical sup-norm estimates established below for wave equations and Klein-Gordon 
equations. 

• Equipped with these refined decay estimates, we are able to improve our initial assumptions 
and close the bootstrap argument. 


^For convenience in the following proof and without loss of generality, we prescribe data at time t = 2. 
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Before we proceed with the details of the proof (which is rather long), the reader may find 
it useful to read through the following heuristic arguments which rely on notations (only briefly 
explained here) to be rigorously introduced only later (in the course of the following three sections). 
Our proof proceeds with a bootstrap argument and considers the largest time interval [2, s*] (in 
the ‘hyperbolic time ’s dehned as = t^ — r"^) within which the following energy bounds hold: 




|J| = fc, |/| + I J| ^ N, 

|/| + |J| ^iV-4, 
\J\ = k, |/| + IJ| ^ N, 
\J\=k, |J| + |J|^Ar-4 


wave / high-order, 
wave / low-order, 
Klein-Gordon / high-order, 
Klein-Gordon / low-order. 


where e,S,Ci are parameters. These bounds concern the energy of the wave component u and 
the Klein-Gordon component v, and distinguish between low-order and high-order derivatives. We 
have denoted by E^n the energy associated with the wave equation (for the flat metric m), while 
are partial derivative operators and are combinations of Lorentz boosts (see below for details). 
The heart of our proof of Theorem 11.11 is proving that, by selecting a sufficiently large constant 
Cl and sufficiently small £,<5 > 0, the above energy bounds in fact imply the following improved 
energy bounds (obtained by replacing Ci by Gi/2): 


E^{s, ^ \cies^\ | J| = fc, 

|/| -f 1 J| ^ TV, 

wave / high-order. 


|/| + |J|^iV-4, 

wave / low-order. 

^ ici£si/2+fe4^ |j| = 

|/| + 1 J| < TV, 

Klein-Gordon / high-order. 

£^,,2 (s, ^ \Cies^Y \J\ = k, 

|/| -h |J| ^ iV-4 

Klein-Gordon / low-order. 


(Of course, it is then a standard matter to deduce from this property that, in fact, s* = -too.) 

To derive the improved energy bounds, we differentiate the equations (|l.lll) with L'’ with 
\I\ + \J\^N-. 

-nd^L-^u = d^L-^ {P^^^davd/sv) + {Rv'^) , 

-nd^L-’v + uH°‘^d^L\ + c^d^L\ = -[d^ ,u H°‘^da,dp]v. 


For these differentiated equations, we perform energy estimates along the hyperboloidal foliation 
and we are led to seek for an integrable time decay for the following the three term^: 

■■=\\d^LJ {P^^d^vdpv) 

(1.14) Ti’\s)-.= \\d^L^ 

For lower-order indices |/| -t | J| ^ N — A, the terms t(’'^{s) and {s) are easily controlled, 
since from the bootstrap assumption and the global Sobolev inequalities on hyperboloids we have 
(basic) decay estimates which lead to time-integrable bounds: 

(1.15) tI\s) + Ti'\s) < s-3/2+(fe+2)5^ provided |/| 4- | J| ^ TV - 4 with | J| = k. 


On the other hand, for higher-order derivatives these basic decay rates are not sufficient and we 
can not conclude directly. In addition, for the third term T^’'^{s) (for arbitrary |J| + |J|), we also 
cannot conclude directly and we need sharper pointwise decay. 

To overcome this challenge, we rely on our sharp decay estimates, established below in 

ProDOsition l3.ll Ifor the wave component) and Proposition [TS] (for the Klein-Gordon component). 


^See II2.7II for the notation. 





PHILIPPE G. LbFLOCH AND YUE MA 


These bounds allow us to improve the basic pointwise estimates, and we find (for all 

\I\ + \J\^N-4): 

\L^u\ < 

Id^L-’v] < C'i£(s/t)2-7^s-3/2+'=‘5, 

\d^L-^dcv\ < 

Returning to our bootstrap assumptions, we thus see that for all |/| + | J| ^ iV 


(1.16) 


^ “ Il+l2 = I 

J2=-J 


< (Ci£)2s-1 + '=^ 

(by assuming, without loss of generality, |/i| + | Ji| ^ fV — 4 in the above calculation). Similarly, 
we also obtain 


(1-17) Wd^L-’ (v^) 11^2(5,^) < (Cie)2s-i+'=^ 

We thus succeed to uniformly control the terms t(’'^{s) and T^''^(s) (for all relevant I, J), and this 
is already sufficient to conclude with the desired improved energy bounds for the wave component. 

Dealing with the last term Tg’“^(s) arising in the equation of the Klein-Gordon component is 
more delicate. Observe that the commutator is a linear combination of the following three types 
of terms: 

{d^^L-’^u)d^^L^^d^dpv, Il+l2 = I, Ji + J2 = J, \Ii \ > 1 , 

(1.18) {L^'^u)d^L^'^dcd0V, J'l + J 2 = J, A > 1, 

udadpd^L^'v, f ^ J-l. 

The first expression above is directly controled thanks to the available sharp decay estimate, while 
for the second and third ones and due to the presence of the term L'^u, a refined decay estimates 
and a Hardy-type inequality (for the hyperboloidal foliation) must be used, as we now explain. 

Let us begin by discussing derivatives of higher-order and consider (for instance) the second type 
of terms in (|1.18|) : for all \ J[\ ^ —4, we use the sharp decay bound combined 

with the energy bound on dadpv (implied by our bootstrap assumption). When |/| -I- 1J 2 I ^ N — A, 
the sharp bound \d^L'^dav\ < Ci£ and Hardy’s inequality are used. We thus 

find 

(1.19) < (Cie)2s-i/2+'=^ 

Dealing with lower-order derivatives is easier and, again, we take the second type of terms in 
(frra as an example: for \ J[\ ^ |/| -f |J| ^ N — A, we apply directly the sharp bound < 

Ciet~^s^^ and the energy bound given by our bootstrap assumption. This leads us to the stronger 
decay 

(1.20) I [Al\ < (Ci£)2s-i+'=^. 

In conclusion, in view of (I1.15I) - (I1.20L we can gain enough time decay for all of the terms arising 
in the evolution of our energy expressions and, therefore, the energy estimate on the hyperboloidal 
foliation leads us to the desired improved energy bounds. 


1.5. A general class of wave-Klein-Gordon systems. The technique presented here applies 
immediately to a much broader class of systems. Indeed, it applies to the following system of 
wave-Klein-Gordon equations 

□ui + B^°‘^UjdadpUi = Fi{u,du,v,dv) = SaUjdpUk + RiV^ + S'^^davdpv, 

□u -I- B^°‘^Ujdadpv — = 0, 

Wi\t=2 = Wio, v\t=2 = Vo, 

^ dtWi\t =2 = Wii, dtv\t =2 = Vi, 


( 1 . 21 ) 
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with unknowns u = (ui) (1 ^ i ^ n) and v defined on Minkowski space while Wig, '^Oi I’l 

are prescribed initial data and c > 0 is a constant. As usual, we assume the symmetry conditions 

( 1 . 22 ) 

and our main assumption is the null condition for the wave components Wi'. 

(1.23) = 0 for all (^ 0 )^ - ^(^a)" = 0. 

a 

In the earlier work [22], the nonlinear terms B^°^^Wjdadpv (actually denoted w^dadpw^ 

therein) were assumed to be vanishing, and in fact this was our only genuine restriction since, 
with such terms, solutions may not have the same time decay and asymptotics of solutions as 
the ones of the homogeneous linear wave-Klein-Gordon equations in Minkowski space. With the 
new technique in the present paper, the hyperboloidal foliation method does extend to encompass 
these terms (provided 5 !“^ is a null quadratic form). 

Let us consider the initial value problem (ll.2ip with sufficiently smooth initial data posed 
on the spacelike hypersurface {t = 2} of constant time and compactly supported in the ball 
{t = 2; |a;| ^ 1}. Under the conditions (ll.22l) - (ll.23l) . there exists a real eg > 0 such that, for all 
initial data WiQ,Wii,VQ,vi : ^ R satisfying the smallness condition 

(1-24) ^ '*^o)IIh'^+ 1(R3) + ll(^^ilj'*^l)llH'^(R3) < COi 

i 

the Cauchy problem (11.211) admits a unique, smooth global-in-time solution. In addition, the 
lower-order energy of the wave components remains globally bounded in time. 

2. The hyperboloidal foliation method 

2.1. The semi-hyperboloidal frame. We begin with basic notions and consider the (3 -I- 1)- 
dimensional Minkowski space with signature (—, -I-, -I-, -I-). In canonical Cartesian coordinates, we 
write {t,x) = and := jeep = {x^Y + addition to the partial 

derivative helds dt = dg and da, we will also use the Lorentz boosts (for a = 1,2,3): 

(2.1) La . X dt tda Xadt) X^da, 

where we raise and lower indices with the Minkowski metric. 

More precisely, throughout, we analyze solutions defined in the interior of the future light cone 

% := {{t, x) / r < t — 1} 

with vertex (1, 0,0, 0), and we introduce the following foliation of the interior of the cone {{t,x) / \x\ < t} 
by hyperboloidal hypersurfaces with hyperbolic radius s: 

■= {(t,x) / Y — Y = s^; t > 0 }. 

The sub-domain of X limited by two hyperboloids (with sg < si) is denoted by 

^bo.sil := {(i: a:) / ^ ^ s?; r < t - 1} ^ X. 

Observe that the hyperboloids eventually “exit” the region X and are asymptote to the same light 
cone {t — r = O}. 

The semi-hyperboloidal frame, as we call it, is defined by rescaling the Lorentz boosts: 

(2.2) do:=dt, d^-=Y^t + da (a = 1,2,3). 

Observe that the vectors generates the tangent space to the hyperboloids. Furthermore, we 
also introduce the vector field d_L := dt -I- ^da, which is orthogonal to the hyperboloids for the 
Minkowski metric. (This vector field also coincides, up to an essential factor 1/t, with the scaling 
vector held S.) 
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To make explicit the change of frame formulas and da 

matrices 


($f) = ($/) 


n 0 

x^/t 1 
x‘^/t 0 
y x^/t 0 


0 o\ 
0 0 
1 0 
0 i) 


K) = 


Hip , we need the following 

/I 0 0 0\ 

-x^/t 1 0 0 

-x‘^/t 0 1 O' 

\-x^/t 0 0 1/ 


Any tensor can be expressed in either the Cartesian natural frame {da} or the semi-hyperboloidal 
frame {da}. We use standard letters for components in the Cartesian frame and we use underlined 
letters for components in the semi-hyperboloidal frame, so that, for example, T“^da ® dp = 
la® dp, and the relations between T“^ and are 

Associated with the semi-hyperboloidal frame, we have the dual semi-hyperboloidal frame 

(2.3) 9° := dt - y t ■= 

and the relations between the semi-hyperboloidal dual frame and the standard dual frame are 
9°‘ = dx’^ = $“/0“ . Hence, for any two-tensor Tapdx'^^dx^ = T_ap^ ®9^, we have the 

change of basis formulas 

Tap = Ta'P'HHp, Tap = Ta,p,HHp. 

With the above notation, in the semi-hyperboloidal frame we can express the Minkowski metric 
and its inverse as 


1—1 

1 

1 

1 

-x^t \ 

—x^jt 

1 — jtY' —x^x’^jt^ 

—x^x^jt^ 

1 

to 

—x^x^jt^ 1 — 

—x^x^ji? , 

CO 

1 

—x^x^jt^ —x^x^jt^ 


/ -s^/e 

—x^jt —x^jt —x^/t \ 


—x^jt 

1 0 o' 


—x^jt 

0 1 0 , 


V 

0 0 1 ) 



Furthermore, given any multi-index I = (a„, q;„_i, ..., ai) (where the order is chosen for 
convenience), we denote by d^ := da^da^-i ■ ■ -^ai the product of n = |/| partial derivatives (with 
0 ^ ^ 3) and, similarly, by L'^ = La^La^_i ■ ■ ■ La^ the product of n = | J| Lorentz boosts (with 

1 ^ Ui ^ 3). 


2.2. The hyperbolic variables and the hyperboloidal frame. Within the future cone DC, we 
introduce the change of variables 

(2.4) = s := — r'^, = a;“, 


together with the corresponding natural frame 


(2.5) 


3 3 '^3 _ 

^0 ^ 

da := 4. = ^dt + da = yd, + d„ 


which we refer to as the hyperboloidal frame. The transition matrices between the hyperboloidal 
frame and the Cartesian frame are 

s/t 0 0 0 

x^/t 1 0 0 

x^jt 0 1 0 I’ 

V x^jt 0 0 1/ 
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= (^f) = = 


( tjs 0 0 0 \ 

-x^js 10 0 

—x^js 0 10 

V -x^js 0 0 1/ 




so that da = ^a^P A = 

The dual hyperboloidal frame then reads dlP := ds = ^dt — ^dx°^ and dx°‘ := dx°‘. The 
Minkowski metric in the hyperboloidal frame read^ 

/ —1 —X^js ——nfii 

—x^js 10 0 


V 


r7s 


-x^js 

0 

1 

0 


-x^/s \ 


/ 


In summary, an arbitrary tensor can be expressed in three different frames: the standard frame 
{da}, the semi-hyperboloidal frame [da], or the hyperboloidal frame {da}- We use symbols, 
underlined symbols, and overlined symbols for tensor components in these frames, respectively. 
For example, a tensor T°‘^da <8> dp is written as 


-a 13 


■a. -r^/3 


where T = 
boloidal frame 

( 2 . 6 ) 


T^^da ®dp = T^^d^®dp = T"^da ® dp, 
and, moreover, by setting C := u\ax.ap \T° 


we have in the hyper- 


-00 


\ ^ c{t/sY 


|r“°|^c(Vs), 


|T I ^ C. 


2.3. Energy estimate on hyperboloids. Throughout this paper, for any function u = u{t,x) 
defined in ]R^+^ (or a subset of it), we consider the integral of on the hyperboloids IKg defined as 
follows: 


(2.7) 


i(„^):=r udx=\ u(Vs2 + r2,a;)da 
^ “ Jar. Jr3 


The subscript refer the fact that we are endowing with the Hat Euclidian metric. We emphasize 
that this is not an integration with respect to the induced Riemannian metric and volume form 
which should be (s/t) dx. This notation will be more convenient for our analysis in this paper but, 
of course, it is completely straightforward to re-state all of our estimates by including the weight 
s/t. 

Consider the hyperboloidal foliation of a region 3C[2,si] = U 2 sSs^si together with the hy- 
perholoidal energy (associated to the Minkowski metric) at some hyperbolic time s e [2, si] 

EmAs, u) ■■= J ({dtuY + ^{dauY + 2{x°-lt)dtudaU + dx 


( 2 . 8 ) 


^ [ ({{sAAtuf + dx 

-■ f ( {dj_uY -f V {{s/t)dauf + V {t~^Vlabu)^ + ) dx, 

\ a a<b J 


where we have also introduced the rotational vector fields flab := x°'db — x^da (not directly used 
here). When c = 0, we also write Em{s,u) := Ejn,o{s,u) for short. 

We will also need the hyperboloidal energy for the curved metric g°‘^ = m“^ -I- h°‘^: 

(2.9) Eg As, u) := E^A^^ J (2h°‘d g^yg^yXa - h°‘^davdpv^ dx, 

where we have used the notation = 1 and = —x°'lt. 

All of our estimates concern the interior of the light cone 3C n {t > 2} away from the origin. 
From here onwards, we assume that all the functions under consideration are spatially compact 


®Our sign convention is opposite to the one in our monograph m, since the metric here has signature 
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and, in particular, vanish identically in a neighborhood of the light cone — 1 = |a;| = r}. More 
precisely, we assume that the initial data on the slice t = 2 are supported in the ball \x\ ^ M for 
some M e (0,1), and we construct solutions supported in the larger domain |a;| ^ M +1. In short, 
we will say that the functions under consideration are spatially compactly supported in X 
or, in short, spatially compactly supported. 

We easily adapt the energy estimate in [22l Proposition 2.3.1] to the equation (11.1111 . as follows. 


Proposition 2.1 (Energy estimate for the hyperboloidal foliation). 1. For every function u which 
is defined in the region X^ 2 ,s] spatially supported in X, one has for all s ^ 2 


( 2 . 10 ) 


E, 


i(s, u)^/^ ^ Em{2, u)^/^ + J 




2. Let V be a solution to the Klein-Gordon equation on a curved spaee (with the definition o/D 
given earlier) 

(2.11) -UgV + c\ = f, 

defined the region X^ 2 ,s] o,nd spatially supported in X. Suppose that h°‘^ = — m“^ satisfies the 

following two conditions (for some constant k ^ 1 and some function M): 

(2.12a) K~'^Eg^c{s,v) ^ Em,c{s,v) ^ K^Eg^c{s,v), 


^M(s)E^,c(s,u)1 /2_ 


(2.12b) J {s/t)(^dah°'^dtvdpv — -dth°‘^davdpv^ dx 

Then, the evolution of the hyperboloidal energy is controlled (for all s ^ 2) by 
(2.13) £;™.c(s,u)1/2 ^ k2£;™(2,u)1/2 + ^2 (11/11^2(5^^.) + M{s))ds. 

Proof. We apply the multiplier dtu to Dit and, by a standard calculation, 

^dt(idtu)^ +2(dau)^) -YjSa{dtudau) = dtuUu. 

a a 

We integrate this identity in X^ 2 ,s] and apply Stokes’ formula by observing that, by assumption, 
the functions under consideration are spatially supported in X, so that there is no “boundary” 
contribution, and we find (see [351 Sec. 2.3]) 

-Em{s,u) --Em{2,u) = I I {s/t)dtu\Judxds. 

^ ^ Ji JMs 

We differentiate this identity with respect to s and apply Cauchy-Schwarz inequality, as follows: 

d C 

Em{s,uy^'^—Em{s,uy^'^ = J {s/t)dtvn\u dx ^ ||□^t||i 2 (J<;^)||(s/^)dt■u||i 2 (J<;^). 

Next, by recalling (12.81) . we find -^Em{s,u)^^‘^ ^ and, by integration over [2,s], the 

inequality (12.1011 is established. 

Next, for the derivation of (12.131) . we rely on the multiplier dtv and, by a standard calculation, 
we get 

^dt{idtvf + ^(dau)^ + c^u^) + ^ da ( - davdtv) 

a a 

+ da( - h'^^dpvdtv) + ]^dt{h°‘^dc,vdpv) 

= dtvf - dag^^^dpvdtv + ^dtg^^^davdpv. 

As was done in the derivation of (12.101) . we combine this identity with (I2.12all - (l2.12bl) and 
obtain 

Eg,c{s,vy^'^-^Eg^c{s,vy^^ ^ niWfh^^^ge,) + M{s))Eg^^{s,vy^^. 
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We integrate this inequality on [2, s] and find 

Eg^c{s,vYl‘^ ^ (II/IIl=(m,)+^( s))ds. 

Finally, we again apply (I2.12a|l . and (I2.13|) is proven. □ 

2.4. Sobolev inequality on hyperboloids. In order to turn an L? energy estimate into an 
L® estimate, we will rely on the following version of the Sobolev inequality (Klainerman |19j . 
Hormander m Lemma 7.6.1], LeFloch and Ma [22l Section 5]). 

Proposition 2.2 (Sobolev-type estimate on hyperboloids). For any sufficiently smooth function 
u = u{s,x) which is defined in 3C[2,+oo) and is spatially supported in %, one has 

(2.14) sup (s + |u(s,a;)| < Ell •)||j;,2(5^^), s > 2, 

(s,a;)6!Xs |/|<2 

where the implied constant is independent of s and u. 

Proof. Recalling that t = s"^ + |a;p on Jfg, we consider the restriction to the hyperboloid 

Ws{x) := u(-\/s‘^ + \x\^,x). 

Fixing So and a point (to.xo) e fifso (with to = -y/sg + |a:op), we observe that 

(2.15) daWsoix) =ffiu{fsl + \x\‘^,x) =ffiu{t,x), 
with t = ^s'o + \x\^ and 

(2.16) tdaWsoix) = tffiu[f s^ + |xp,t) = Lau{t,x). 

We introduce the function gso,to{y) ■= Wsg(xo + tg y) and write 

5so.to(0) = Wsoixo) = u{f sg 4- |a;oP,xo) = u{to,xo). 

By applying the standard Sobolev inequality to the function Qso.to, we find 

|ff.o.to(0)|'<c 2 r \d^9soM?dy. 

|/|^2 

where B(0,1/3) c is the ball centered at the origin and with radius 1/3. 

Next, taking into account the identity (with x = Xg + tgy) 

da9so,to{y) = todaWsoixg + tgy) 

= igdaWsoffi) = tgffiu{t,x) 

and, in view of (I2.15|) . we find (for all I) dffiso,taiy) = {tod)^u{t,x) and, thus, 

|5^o,to(0)|^ E f \itoiyu{t,x))\^dy 
|/|^2 

=Cto^ E f \{todYu{t,x))fdx. 

|7|sg2 •lB{{to,xo),to/3)n^sQ 

We note that 

(tgffiitoffiWso)) =tlffidf,Wso 

= {to/ty{tffi){tffi)wso - {to/tf{x'^/t)Li,Ws„ 

and that x°‘/t = Xo/t + ytg/t = (ccg/to + y){to/t)- Consequently, in the region ye R (0,1/3) of 
interest, the factor |a;“/t| is bounded by C{to/t) and we conclude that (for |/| ^ 2) 

litgdYul sS 2 \L'^u\{tg/tY. 
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On the other hand, in the region jaiol < io/2, we have to ^ and thus 

to ^ Cso ^ C'\J\x\^ + Sq = Ct 

for some fixed constant C > 0. When |xo| > to/2 then in the region i3((to,xo),to/3) n Jfso we 
have to ^ Clx] ^ Cy/\x\‘^ + Sq = Ct and, consequently, 


Kto^^ul^C 2 \L^u\ 


and 


5so.io(yo)|^ ^C'to^ ^ \{tdYu{t,x))'\ dx 

|/|sS 2 •t-B(xo,to/ 3 )nXCsg 

^Ctg^ Xj f \L^u{t,x)\^ dx. 

|/|!S2 ‘^^'>0 


□ 


2.5. Hardy-type estimate along the hyperboloidal foliation. The following variant of Hardy’s 
inequality was established in [22l Section 5]. This inequality plays an essential role in order to 
estimate the norm of the wave component itself (but not only its gradient). 


Proposition 2.3 (Hardy-type estimate on the hyperboloidal foliation). For any sufficiently 
smooth function which is defined in the future region ^[ 2 ,s] o,nd is spatially supported in OC, one 
has for s > 2 

Is ^ \\u\\l^{K2) + Xj 

(2.17) rs _ J 

+ X ^ \\i^/Y^aU\\Lj{Ks)) ds, 

a ‘.'2 

where the implied constant is independent of s and u. 


The proof uses a version of the classical Hardy inequality on hyperboloids, as well as a vector 
field that will be introduced in the proof of the proposition, below. 


Lemma 2.4. For any sufficiently smooth function which is defined in the future region ^[ 2 ,s] 
is spatially supported in 3C, one has for all s ^ 2 

Ik ~ X 

a 

where the implied constant is independent of s and u. 


Proof. As in the proof of Proposition 12.21 we consider the function Ws{x) := u('y/ s^ + |xp, x ), 
which satisfies daWs{x) = + kp, x), and we apply the classical Hardy inequality to Wg. 

It follows that 


r \r ^Ws{x)\'^dx < r |Vws(x)kdx = CX f |^a“(V+ r^,x)|^dx 
Jb3 Jb3 „ Jb3 


< 


sL 




□ 


Proof of Proposition \2.S\ Let y be a a smooth cut-off function satisfying 


X{r) 


0, 0 sS r s: 1/3 

1, 2/3 ^ r. 


and let us distinguish between the region “near” and “away” from the light cone. We consider the 
decomposition 
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Our estimate of ||(1 — is based on the inequality (l — x(r/t))s ^ ^ Ct ^ so 

that, by Lemma [2.41 


(2.18) 


1(1 - X(r/f))us ^||L2(jr,) < l|i 


s;Ik < cY,\\s.a 


u r2 




The estimate near the light cone is more delicate and we first observe that, in the region 3C[2,, 
of interest, x{r/t) < and, thus. 


\\x{r/t)i 


< C'lk(l +'’^) ^^^X{r/t)s ^u||i2^), 


and the right-hand side of this inequality is controlled as follows. We introduce the vector field 
W = (O, —and compute its divergence 


diyW = -2s~^d a;“tx(r/t) rx{r/t)u x'{r/t)r 


(2.19) 


(1 -I- r'^yl'^s r(l 4- 


2s- 


r s(l + (1 -I- 7 ’ 2 ^i /2 


r^t 4- 3< 


4- 2- 




,(l4-r2)2s2 (l4-r2)s4 

By applying Stokes’ theorem in the region dC[ 2 ,si]j we find 


wxkA))^- 


I 


W ■ nda 


div W dxdt =1 IT • n dcr 4- [ 

■JD<s 

f Y^\^x{r/t)s-^\'dx - { 

i -I- r jj 


212 1 + 


Differentiating this identity with respect to s leads us to 


( 2 . 20 ) 


d/r 

ds V Ja 


div IT dxdt ) = 


-ff 

ds V 


r I 
1 4- j,2 I 


ix{r/t)i 


ix{r/t)i 


dx 


dx. 


= 2 


ruxirjt) 

d 

rux{r/t) 

s(l 4- r2)i/2 

L2(TC,) ds 

s(l + r2)i/2 


L^i^s) 


We then integrate (12.191) in the region 3C[2,si] c IK n {2 < ^/t'^ — r'^ ^ si}: 

rx{r/t)u x‘^tx{rlt) 


I 


div IT dxdt = — 2 


s 

^[2,31] 


-1 


daU- 


(l4-r2)i/2s^(l + r2)i/2 

i'U rx{rlt)u x'{rlt)r 


dxdt 


i 


2C[2,.P rs(l + r2)l/2(l + ^2)l/2 


dxdt 


r^t + 3t 


+ 2 - 


r^t 


3C[2,.i]V(H-( l4-r2)s4 


[uxir/t)Ydxdt, 


which yields the identity 


I 

•^3C[2.n] 


div IT dxdt = — 2 


p r {s/ty 
J2 J:k,, 

rsi 


-1 


daU 


rx{r/t)u xHxir/t) 

(l4-r2)i/2sr(l-4r2)i/2 


dxds 


dxds 


r / /.N - 1 ^ rx{rlt)u x'{r/t)r 

J2 r s(l 4-r2) 1/2 (1 + 7-2 )i/2 

r t i\f r^t + 3t ^ r^t , ^^2 , , 

( n ^.2W,2 + 2 ,1 ^.2 W4 ) i^xir/t)) dxds 

rsi 


(1 4- r2)2s2 (2 4- r2)s‘* 

= : J ' (Ti + T2 -f T3) ds. 
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Here, we have T 3 ^ 0, while 


Ti = - 2 s- 


^ 2 s" 


ruxirjt) 


'(l + r2)i/2sr(l + r2)i/2 




s(l + r 2 )i /2 




^ ||(s/t)daW||L2(Jt,)||x(?’A)a;“t?’ ^(l + r^) 1/2|| 




and 


^Cs- 


r, = - 2s- 


^Cs- 


^Cs" 


iXirlt) 


^||(s/i)daw||i 2 ^) 


s(l + r2)l/2 

M rx{r/t)u x'{r/t)r 


f (^A)- 


Jm, 'r s(l + r2)i/2 (l + r2)i/2 

rux{r/t) 


dx 


s(l + r2)i/2 
ruxirjt) 




ur ^\\L^Mj\rx'{r/t){l + r^) 




a 




s(l + r 2 )i /2 

where Lemma 12.41 was used. 

We write our identity in the form ^ ^ div W dxdt^ = Ti + T 2 + T 3 and obtain 


( 2 . 21 ) 


_d 

ds 

^ s 


-1 


div W dxdt 

‘^[ 2 , 31 ] 

ruxir/t) 


s(l + r2)i/2 


a 


,) -t- ll^a“llL2(M^) 




Finally, combining (12.2111 and (12.2011 yields us 


_d 

ds 


^X{rlt) 


s(l + r 2 )i /2 
and, by integration over [ 2 , s], 


L 2 (aCs) 




( 2 . 22 ) 


||r(l + r2) i/2^(r/t)s 

\\l^j<2 ) ' Zj 


sg ||r(l + r^)“^/^x(r/t)2“^u|L„..^ , + Y! f s“Vll-daU 




From Lemma [2.41 we then deduce that 

||x(r/t)s"^M||i 2 (;jc^) < \\r{l + r^)-^^^x{r/t)s-^u\\L 2 ^^^) 

(2.23) 


^ l| 2 "^w 


L 2 (M 2 ) 


+ Xj + ll^aW||L2(j^_))ds, 


and remains to combine (12.1811 with (12.2311 . 


□ 


3. Sup-norm estimates for the wave and Klein-Gordon equations 

3.1. Statement for the wave equation. 

Proposition 3.1 (A sup-norm estimate for the wave equation with source). Let u be a spatially 
compactly supported to the wave equation 

-□« = /, 

u\t=2 = 0 , dtu\t=2 = 0 , 


(3.1) 
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where the source f is spatially compactly supported in X and satisfies the estimate 

\f\^Cft-^-^t-r)-^+^ 

for some constants Cf > 0, 0 < fi ^ 1/2, and 0 < |i/| ^ 1/2. Then, the following estimate holds: 

i^it-rY-''t-^, 0<v^l/2, 

(3.2) \uit,x)\ < / "c/w ^ 1 , 

\^(t-rrt-^-\ -l/2^iy<0. 

We recall that the energy estimate on wave equation does not control the solution itself but 
only its gradient. So when we apply the Sobolev inequality and obtain a sup-norm estimate (cf. for 
example [22] ), there is no immediate estimate on the sup-norm of the solution itself. The estimate 
above yields a (sharp) sup-norm estimate on the solution itself and will play an essential role for 
the control of the quasi-linear term udtdtv in our model problem. We emphasize that the range 
—1/2 ^ 1 / < 0 will only be used in the second part |23) . 

3.2. Proof of the sup-norm estimate for the wave equation. We now state a technical 
lemma and give the proof of Proposition 13.11 but postpone the proof of the lemma to the end of 
this section. Let da be the Lebesgue measure on the sphere {|?/| = 1 — A} and x e with r = |x|. 
We are interested in controling the integral 

da{y) 


J(A) = I{X,t,x/t) := 




%| = l-A,|f-y|sSA-t 1 {^~\j~y\) 

Our bounds below are consistent with the obvious estimate where x = 0: 
(3.3) /(A, t, 0) = 47r(2A - l)-i+'^(l - A)^. 

Clearly, when 0 < A ^ /(A) = 0. 

Lemma 3.2. When ^ A ^ 1, the following estimate holds: 


i-fj. ■ 


IW ^ 


At(l — A) f t — r 

p,r 

(1-A) 


t + r 


- A 2A 


t + r 


—1 + /2 


t — r + 1 , t + r + 1 

t + r + 1 , t — r 

< A ^ , 


2t 


. , , t -l- r -t- 1 t ■ 
provided -^ 


(1 — X)t ft — r 
p.r V ^ 

Proof of Proposition From the explicit expression 


max 


t — r t + r + 1 


2t 


2t 

sg A ^ 1. 


(3.4) 


= ^ \ f f{s,x-y)dads 

477 J2 t S 


in which the integration is made on the intersection of the cone {{s,y) / \y — x\ = t — s,2 ^ s ^ t} 


and {(t,x) /r < t— \,t^ 


s^, t > 2 }, we obtain 


|u(t,x)| ^ 


477 


n 

J2 J|' 


s ^ ’'{s — |x — y\) 


dads 


2 

Cf 


|y|=t-s,|a:-y|sgs-l 
1 


477^1 + '^-'' 

Cf 


rf 

Jf J|y'| = i-A.|f-j/'|s:A-t- 

p(i-A)-iA-2-‘^ r 

h hv'\- 


t — s 

(1 - X)-^X-‘^-''dadX 


(A. 


y'\) 


l-M 


(A := s/t, y' := y/t) 

dX. 


^ {'^ —If — 2 /^ 1 ) ** 

When \ j — y'\ P X — t~^ holds, we obtain ^ A sS 1. For convenience in the notation, in the 

following calculation we replace y' by y. We first assume that r > 0 and we distinguish between 
two main cases: 
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Case 1: r ^ In Lemma 15^ all three cases are possible: 


\u{t,x)\ §— r (1 - a)-ia- 2-" r 


da 


|y| = l-A,|f-y|s:A-t-i (A - |f -y|) 


1-M 


dX 


<- 


Cf 




(l-A)-iA 


_i,_2-^At(l - A) f t-^ 


dX 




— l+/i 


+ 


Ct 




r (i-A)-i 

Ji=i. 


-1,-2- - A)t f t-r 


dX, 


thus 


\u{t,x)\ < 


Cf t ft 


r V t 


x-^-'dx 


+ 


£L- p" ^-2-uft±l 

Jl 


X- 


C + ^-P- L + r + 1 I, t 


-A 2A- 


t + r 


— 1+l-L 


dX 


+ 


Cf t f t — r 


A* /-l 


x-^-'' dX. 


fiC+''-^r\ t ) J^r. 

For the first integral, we recall that r ^ and that 0 < |ia| ^ 1/2, and write 


t + r + l 


t —r+1 


X-^-''dX < 


t — r 


i+i' 


< 1 , 


so that 


Cf t ft — r 


r t 


A"^"^ dX 


< Cffi~\t- r)>^t 


- 1 - 1 / 


Next, for the second integral in the right-hand-side, we just remark that 


t 

L 


\ —2 — 1 / 


t + r 


-A 2A - 


t + r 


— 1 + /^ 


< 


t — 

Jill 


t + r+1 


2A- 


t + r 


— 1 + /J. 


dA = - 2A 


dX 


t + r 


~ 1 

< —. 

t + r + l fL 


This leads to 


Cf 

tl + u- 


J< 


A 


t + r + 1 


-ir 


-A 2A - 


t + r 


-i-ffj. 


dX < 


Cf 




For the third term, in view of ^ ^ 1, we obtain 


Cf t ft 




x-'^-'dx <- 


Cf t ft- 


M rl 


dX 


dX 


j. y f J ^fi+v-fi j, y f j 1 1 -,. 

<Cffi-\t-r)>^t-^-'. 

So we conclude that in the case 0 < r ^ \u{t,x)\ < Cffi~^{t — r)^^t~^~''. 

Case 2: ^ o r ^ The second case in Lemma 13.21 is not possible, and we have 


\u{t,x)\ < 


Cf ft 


I t 


r 

I t + r + l 

2t 


X-^-’' dX + dX . 
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Since ^ second integral is bounded by a constant C. For the first integral, we see 

that when i/ > 0, 


'd\ < 


1 f t — r + 1 


— V 


So in this case when z/ > 0, we obtain \u{t,x)\ < Cf{^xv) — r)^ ^. 

When V < 0, we write 




dA < —r 
\v\ 


1 / t + r + 1 


— V 


1 



and, therefore, we obtain |M(t,x)| < C/(/i|v|)“^(t — . 

When r = 0, we make the following direct calculation, remark that in this case, ^ A sS 1, 
by (l33t : 


\u{t, a;)| ^ ^ - r (1 — A) ^A ^ dA r 


da 


< 


Jl - N-U-»-(2 - A)-‘+<'(l - A)VA 

(since r = 0), 


M A* 


which completes the proof. 


□ 


Proof of Lemma \S.2\ When r = 0, the estimate is trivial. When r > 0, without loss of generality, 
let X = (r, 0,0). The surface S\ := {|y| = 1 — A} n {|| — y\ ^ A — t~Y is parameterized as follows: 

• 9: angle from (1,0, 0) to y with 0 sS d ^ tt, 

• (j): angle from the plane determined by (1, 0,0) and (0,1,0) to the plane determined by y 

and (1, 0,0), with 0 ^ ^ 27r. 

Then, we have y = {I — A)(cosd,sindcos(/>, sindsin;/)) and we distinguish between two cases, as 
follows. 


Case 1. When * ^ A ^ , we only have a part of the sphere {|y| = 1 — A} con- 


2t 


2t ’ 
f-ll 


tained in the ball {|f — < A — t where cos(d) > 

'(r/i)" + (l-A)=-(A-t-i)' 

(2r/t)(l-A) 


(r/t) = + (l-A) = -(A-t-i)' 
(2r/t)(l-A) 


. So we set do := 


arccos 


and see that 


- y = A- 


— + (1 — A)2 — 2- (1 — A) cos d 


and dcr = (1 — A)^ sin(d)ddd^. The integral is estimated as follows: 


da 


iiy| = l-A.|f-y|sgA-t-i (A— If — j/|) ^ 

= J d(j)^ (1 — A)^ sind^A — \J'^ + Y — A)^ — 2^(1 — A) cosd^ 


— l + /i 


dd 


— 1+/J, 


=27rJ (1 — XY sind^A — + (1 — A)2 — 2^(1 — A) cos d^ 

r^Q / j 2 \ + 

= — 27r(l — A)^ J ^A — Y ^ + (1 “ A)2 — 2^(1 — A) cos dj dcosd 
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thus by setting oj = cos 9 


r da 

J|y|=i-A.lf-y|=sA-t-i (A - If - y\f~^ 


V 2 \ 

^ + (1-A)2-2^(1-A)a;j dec 


7rt(l — A) 


_i+. ,^t(i-A) 




(A-V 7 ) = 2 


r'+'^(A-C)dC, 


where we have used 7 = f 2 - + (l — A)^ — 2f(l — A)a; and C : = A — Then, we distinguish between 
the following two sub-cases. 

Case 1.1: f ^ 1 — A or, equivalently, A ^ We now find 


,7rt(l-A) r^-l7-(i-^)l 


r .f-i 


ri+'^(A-C)dC 


rt(l-A) 


C^+>^{\-QdC,< 


Case 1.2: 1 — A < f or, equivalently, A > We find 


At(l — A) (< — rY 

fir 


rt(l-A) r^-l7-(i-^)l 


C-'+^(A-C)dC 


^7Tto_^ r ‘ ^_i+^ ^ 

r Jt-i 


At(l - A) (f - rY 

fir 


Case 2. When ^ A ^ 1, the sphere {|i/| = 1 — A} is entirely contained in {\{x/t) — y| ^ 

A — t~^}: 

r da C da 


r da _ 1 

r 

J|y| = l-A.|f-j/|=SA-t-i (A- 

f-yy-^ J 

1 

1 

7\ 


271J (1 —A)^sin0^A 


I 2 \ —1 + /^ 

(l-A)2-2^(l-A)cos6»j de 


= 241 -A)M fA-y^ + (l-A)2-2^(l-A)wj dec 


and thus 


r da 2 ^^^^ ~ 

J|!/| = l-A,|f-y|s:A-t-i (A-If r 


,A-If-(1-A)| 


r J2A- 


r'+^(A-c)rfc 

'■ JA-(f + (l-A)) 

TTtfl - Al r^“l7-(l-^)l 

^ ^ C^+^{X-0dC. 

J2\-i±X- 


We now distinguish between two sub-cases. 


Case 2.1: When f ^ 1 — A or, equivalently, A ^ ^-r^, we find 


, 7 rt(l-A) r^-lf-(i-^)l 


r'+^A-odc 


rt(l-A) 


r'+^(A-C)dC<C(l-A) 


-A 2A - 
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where we have observed that in the integral the function is decreasing and we can 

bound this integral by the value of the function taken at the inferior boundary (which is 2A— 
times the length of the interval which is 2 r/t. 

Case 2.2: When 1 — A < 7 or, equivalently, A > we have 


A-|r_(i_A)| 


C-'+''(A-C)dC 


^ 7rt(l - A) r 
r h\-^^ 

r'WA-OAc 

r J2\-i±^ 


t — r 

c'(i-A)- r ‘ c^+f^dc 

r J2\-1 ±l 




fir 


C{l-X)t (t-r 


fir 


t+r+l ^ t-r 
2t ^ t ’ 

inequality holds true. 


When ‘+^+1 ^ both case above may occur, while only Case 2.2 is possible if the opposite 

□ 


3.3. Statement for the Klein-Gordon equation. Consider (sufficiently smooth and spatially 
compactly supported) solutions to a Klein-Gordon equation on a curved space and, specifically. 


(3.5) 


- Clgn -f c^v = /, 

v\k 2 = Vo, dtv\ji^=vi, 


with initial data Vo,Vi given on IK 2 and compactly supported in ^2 and the metric has the 
form g°‘^ = with is spatially compactly supported in % with sup \ h | ^ 1/3. 

Before we can state our estimate, we need some notation. Given a constant C > 0 and using 
the notation s = — r^, we consider the function 

ht,xiX) := /i°'’(At/s, Ax/s), 

and, by denoting by for the derivative with respect to A, 

K xW = -dth^^{\t/s,Xx/s) + —dah^'^{Xt/s,Xx/s) 

’ s s 

= -d_i_ ft,™ (At/s, Ax/s). 


We set 


(3.6) 


So := 


t -I- r 
t — r' 


O^r/m 3/5, 
3/5 ^ r/t ^ 1, 


and introduce the following function V which is defined by distinguishing between the regions 
“near” and “far” from the light cone: 


K:= / 


(lko||L»(M2) + I|i’iIIl”(M2))(i + 

-fi^(s)-f J F(s)|/i/^(A)|e'^S=IG.»G)|dA^- q sS r/t ^ 3/5, 
F(s) H- r F(s)|ti;_,,(s)|e'^S.G?A/JA)|dA^- 3/5 ^ ^ 

Jsn 


F{s) := J ^(i?i[x] + i? 2 [r] + i? 3 [r])(At/s, Ax/s)-I-A^/^/(At/s, Ax/s)^(tA 


with 
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and 


;= s^/'^Y^dadaV + 

a 


~s^ 


dadbV 


+ 


3 v-t 


R2[v\ ■■= + 3s^/^(5oi;^ - {2li^ d^dbV + h^^dadbV + Sp^v), 

i?3[p] := 7i°° ^2x“s^/^dodo'i; + ‘^^'SaV + . 

With these notations, our result is as follows. 


Proposition 3.3 (A sup-norm estimate for the Klein-Gordon equation with source). Considering 
the Klein-Gordon problem (1^ for a every sufficiently smooth and spatially compactly supported 
solution V defined the future region 3C[2.+oo), one has (for all relevant {t,x)) 

(3.7) a;)| -I- {s/t)~^s^^'^\dj_v{t,x)\ < V{t,x), 

This result is motivated by a pioneering work by Klainerman m and the decomposition in 
Lemma l3.4l below. An analogue statement in two spatial dimensions and flat Minkowski spacetime 
is discussed in [28]; see also the earlier work [9]. 

3.4. Proof of the sup-norm estimate for the Klein-Gordon equation. We begin with two 
technical results. 


Lemma 3.4 (A decomposition identity). For every sufficiently smooth solution v to (13.511 . the 
function 


wt,x{^) ■= A^/^u(At/s, Ax/s), {t, x) e X, 
satisfies the following second-order ODE in A 


+ - =00 - 

1 + /i (At/s, Ax/s) 




= (l + h°°(At/s, Ax/s)) ^(i?i[u] -f i? 2 [u] + Ralv] -\- s^^^/)(At/s, Ax/s). 


Lemma 3.5 (Technical ODE estimate). Let G be a function defined on an interval [so,si] and 
satisfying sup|G| ^ 1/3. and k be an integrable function defined on [so,si]. Then, the solution z 
to the ordinary differential equation 


(3.8) 


-z(A) = fc(A), 


+ —C(A)' 

2;(so) = Zq, 2 :'(so) = zi, 


with prescribed initial data zo,zi satisfies the uniform bound 

(3.9) |z(s)| + |z'(s)| < (|zo| + ki|+df(s))+ f (|zo| + |^i| + i^(s))|G'(s)|e^KI«'(^)l‘^^ds 

Jsq ^ ^ 

for all s e [sq, si] and with K{s) := |fc(s)| ds and a suitable constant G > 0. 


Proof of Lemma \3f4[ 1- Decomposition of the flat wave operator. By recalling s = Vt^ — 
and r = |x|, an elementary calculation shows that the flat wave operator □ in the hyperboloidal 
frame reads 

(3.10) - □ = dodo - ^dada + 2^ ^doda + ^do. 

a a 

Given a function v, we can set 

w(t,x) = s^/^w(t,x) = A ~ \x\‘^)^^‘^'v{t,x), 


and obtain 
(3.11) 


s^/^Dw = dodow - ^dadaW -f 2^ ^dodaW - - Xj 
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Again, we define a function of a single variable by 

wt,xW ■= w{Xt/s, Xx/s) = X^^^v{Xt/s, Xx/s). 

We see that 


-^Wt,x{X) = (do + s ^x°‘da)w{Xt/s,Xx/s) = -dj_'w{Xt/s,Xx/s) 
ctA s 


and 

(3.12) 


X ^ ^ X X -:r ^ 


-^Wt,x{X) = ( dodo + 2—doda + -^dadb ]w{Xt/s,Xx/s). 


Combining this with (13.lip and recalling that w{t,x) = s^'^i;(t,x), we get 


(3.13) 


dodo + 2 —dod, + ^dadfc ]w 


x°^x^ 


= + ^dadaW + ^^-^dadbW + 2 -^^aW = 


2. Decomposition of the curved wave operator. We write 

—□« = h°‘^dadpv — (?v + / 
and, by performing a change of frame, 

h°‘^d^dpV dp'V 

=/i°°dodon + 2h^^dodbV + h'^^dadbv + /i“^dad>^ dp'V. 

Then, we obtain 

—= — s^^^h^^dodov — dodbV + h°‘'^dadbV + h°‘^da'^p dp'v) — 


= — h'^^dodo[s^^^v) — 
+ 


( 47 ^ 3 s^/^don^ - dodbV + h°'^dadbV + dp^v) + 


and we conclude with 

«3/2|— I., ^ 


— S 


□a = —h^^dodow — c^w + 3s^'^^don^ 

— s^/^( 2 /i°^dod{,n + h°'^dadbV + h°‘^da'^p dp'v) + 
= —/i°°dodow — c^w + R2 [v] + 

We then combine (|3.13l) with (13.141) and obtain 


(3.15) 




x°'x^^ -^00^ ^ 


dodow + 2—dodaW + T dadbw — h""dodow + c^w = i?i[u] + i?2[n] + 


3. Conclusion. We continue with p.lSp and write 


(1 + h°°) ( dodo + 2—doda + ^^-^dadb )w + c' 


XX 

~ 


and, so, we have 


T-oo/^a;®^^ x'^x^ 

— h I 2—doda H-5— 


X*^ x^ x^ 

dodo + 2—doda H- yp—dadb jw + 

tOOn-1 


dadb]w + i?i[u] + i?2[n] + 


1 + h 


■00 


= {l + h ) \R,[v] + R2[v] + Rs[v] + s^^^f). 


(3.16) 
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It follows that 


(3.17) 


(P „„ , Pwt,x{X) 

1 h {Xt/s, Xx/s) 

= [l + h°^{Xt/s,Xx/s)) ^{Ri[v] + R 2 [v] + Rslv] + f){Xt/s, Xx/s). 


□ 


Proof of Lemma \3.5\ We simply need to integrate out the ODE. We consider the vector field 
h{X) = ( 2 ;(A), z'(A))^ and the matrix A(A) := ^ ^ G)~^ 0 ) ' b' = Ab+ ^ ^ ^ • 

and introduce the diagonalization of A = PQP~^ with 


and 



Q = 

/ m(l + G) 

0 


V 0 

-ic{l + G) 


f 1 

1 \ 

, ( 

p = 

1 ic 

ic 1 5 

P -1 = 


V (1 + G)1/2 

(l+G)l/2 J 

V 


-1/2 


1/2 \ 

1/2 ) 

We thus have b' = PQP~^b + ^ ^ ^ > which leads us to 

{p-^)' = Q{p-^b) + {p-^yb+p-^(^iy 

We regard the term [P~^)'b as a source term and, by a standard formula, 

p-i&(s) = + J eK Qis)dsp-i ° ^ dA 

+ r eKQ(«d'^s(p-i)'(A)6(A)dA. 

'JSq 

Recall that when sup;^£|-i_g] |G(A)| ^ 1/3, the norm of P{X) and P~^{X) are bounded for A e [sq, s]- 

We also remarks that the norm of (P”^)^(A) is bounded by C'|G"(A)| with C a constant depending 
only on c, and the norm of Q is also bounded by a constant C > 0. Furthermore, we observe that 


J Q{s)ds = 


5 _ / *cJ^(l + G) ^/‘^{s)ds 


0 


and thus 


SJ Q(s)ds ^ 


ei\ 


0 —fcj^(l + G) ^^'^{s)ds 

gicSJ(l + G)“i/^(s)ds Q A 


^ 0 g-icSJ(l+G)-i/^(s)ds j ■ 

The norm of the matrix Qi>^)ds jg uniformly bounded by a constant, and the following estimate 
is now proven: 

|z(s)| + |z'(s)| ^ G(|z(so)| + W{so)\) + CK{s) + C f |G'(A)|(|z(A)| + |z'(A)|) dX, 


and we conclude with Gronwall’s lemma. 


□ 


Proof of Proposition \S.A The proof is based on a combination of the bounds (13.91) and (13.171) . By 
recalling the definition of Wt,x{X), we have 

Wt,xiX) = X^I'^v{Xt/s,Xx/s), 

w/^{X) = ^A^/^u(At/s, Ax/s) + -A^/^d_i_u(At/s, Ax/s). 
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That is, Wt^x is the restriction of a:) = s^/'^v{t,x) on the line segment {(At/s, Ax/s), A e [so,s]}. 
We then apply p.9l) and (I3.17P to this line segment, with 

2, 0 sS r/t ^ 3/5, 


So = 


jt + r 
t — r' 


3/5 ^ r/t ^ 1. 


This segment is the part of the line {(At/s, Ax/s)} between the point (t,x) and the boundary of 

-^[so, + oo)- 

Recall that v is supported in %. and the restriction of v on the initial hyperboloid IK 2 is 
supported in 5f2 n X. We recall that when 3/5 ^ r/t ^ 1, wt,x{so) = 0 and when 0 ^ r/t ^ 3/5, 
u’t,x(so) is determined by vg. 

When 0 ^ r/t ^ 3/5, we apply dSSl) with So = 2. When A = 2, we write Wt,x(2) = 
w(2t/s,2x/s) = 2^/^n(2t/s, 2x/s) = 2^/^ro(2x/s), and 

d_ 

3V2 


^'s,x(2) =^{>'^^^v{Xt/s, Ax/s)) 

'-v{2t/s, 2x/s) + 2^/^(s/t)“^dxu(2t/s, 2x/s) 

~-v{2t/s,2x/s) + 2^/^(s/t)“^(5tu(2t/s, 2x/s) + 2^/^(x“/s)dar(2t/s, 2x/s) 
-Uo(2x/s) + 2^/‘^{x°' / s)daVg{2x/s) + 2^/^(s/t) ^Ui(2t/s, 2x/s). 


2 

3 V 2 

3 V 2 


Recall that when 0 sS r/t ^ 3/5, we have 4/5 ^ s/t ^ 1. So we see that |wt,a;(so)| + \w[ 3 ,(so)| ^ 
C(ll'yollL«>(M 2 ) + lkillL«>(M 2 ))- Then by and (I3.17P we have 

|wt, 2 ;(s)| + |w/ 2 )(s)| < C'(||^'o||L®(M 2 ) + l|^'l||L<»(S^ 2 )) + CF{s) 

+ C'(II^^o||l'»(3C2) + II^^i|Il‘®(M2)) J \K,xis)\^^ 


+ c 


£F(s-)|hL(s-)| 


eCV,KxW\d>~ ds. 


We recall that when 3/5 ^ r/t ^ 1, Wt,x(so) = dv/^(so) = 0 and so we have 

kt,.(s)| + kCWI <^^^^(5) + C r ns)|ti;,,(s)|e^«l'*LWM^ds, 

Jsq 

which leads to |ryt,x(s)| + |rc( 3 ,(s)| < V(t,x). It remains to recall the relation between v and w, 
that is, v(t,x) = s^/'^wt,x{s) and 

{s/t)-^s^/'^d^v{t,x) = w/^{s) - ‘^s^/^v{t,x) = w/^{s) - ^s-^wt,xis), 


and the desired estimate is established. 


□ 


4. Commutator estimates 

4.1. Algebraic decomposition of the commutators. We consider the commutators [X, F]w := 
X{Yu) — Y{Xu) of operators associated with our vector fields when the function u is defined in 
the future cone X = {|x| < t — 1}. Our uniform bounds rely on homogeneity arguments and on 
the observation that the coefficients of our decompositions are smooth in X. 

First of all, the vector fields da, and La are Killing fields for the (flat) wave operator □, so that 
the following commutation relations hold: 

(4.1) [da, □] = 0, [La, □] = 0. 

By introducing the notation 

(4.2) [La,dp\=-.Qlpd^, [da,dp\=-.t-^Ml^d^, [La,dp\=-.^f,d^, 
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Mlo = 0, 




Ml, = 


we find easily that 

©lo = -<^2, 

(4.3) Ml, = -^^5Z = ^15Z^ 

02o = -<52 + y<5j = -<52 + $S<5o", 

where $ and 4> were dehned at the beginning of Section [H All of these coefficients are smooth in 
the (open) cone X and homogeneous of degree 0. Furthermore, we can also check that 

(4.4) 0a& = O, so that [La,df,] = 

which means that the commutator of a “good” derivative dj with La is again a “good” derivative. 
(That is, these derivatives enjoy better decay compared to the gradient itself.) 

Lemma 4.1 (Algebraic decomposition of commutators. I). There exist constants such that 

( 4 . 5 ) [dlLa]= 2 

Proof. We proceed by induction on |/|. For |/| = 1, the result is guaranteed by (14.21) . Suppose 
that (1431) holds for all |/i| ^ m, we will prove that it is still valid for |/| ^ m + 1. Let / = 
{a, am, Cim-i, ■ ■ ■ ,cti) and Ii = (am, C(m-i, ■ ■ ■, OLi), so that . Then we have 

[dlLa] = = d„([d^bLj) + =dj Y, - ©L^7<5"^ 

= I] 

\J\^\L\ 

which yields the desired statement for |/| = m + 1. □ 

Lemma 4.2 (Algebraic decomposition of commutators. II). There exist constants such that 

(4.6) [Y,dJ= O^d^Ll 

Proof. We proceed by induction and observe that the case |/| = 1 is already covered by (14.21) . We 
assume that (1^ is valid for |/| ^ m and we will prove that it is still valid when |/| = m+ 1. For 
this purpose, we take = LaL^^ with |/i| = m, and we have 

[Lida.] =[LaL^\da.] = La{[M\da.]) + 

= La( Y +E©L<57Y^ 

^ |j|^|/i|-i,7 ^ 'y 

= E OilLad.L-’ + 

MI=Sbl|-1.7 7 

and, therefore, 

[Llda.]= Y 0l]d^LaM+ Y 0il[La,d^]M + YQZo.3M'^ 

|./|s:Fi|-1,7 |J|sS|/i|-1.7 7 

= E 0lld,LaJ'’+ Y ^«7©i;<5yL'^+E©L<57YL 

\J\^\Il\-ln |J|=S|/i|-1.7 7 

□ 


An immediate consequence of (14.61) is 

(4.7) [d^Llda.]u= E oYj.d^dM-’'u. 

\.J'\<\-na 
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Lemma 4.3 (Algebraic decomposition of commutators. III). In the future cone OC, the following 
identity holds: 

(4.8) [d^L\d^]= 2 L^', 

|R|sO|,|/'|S|/| 

|J'| + |J'l<|r| + |J| 


where the coefficients are smooth functions and satisfy (in 3C) 


(4.9) 




C{\Il\J\,\h\,\Ji\)t-\^^\ 

cd/UJU/iUJiDH^^i-i 


when \J'\ < |J|, 
when |/'| < |/|. 


Proof. Consider the identity 


|/ll + |Jll<|/| + |J| 


In the first sum, we commute and dj and obtain 


[d^pJ^dp] =^}[d^L\d^] 

+ 2 dPLP^JdjdPpP + 2 dPLP^}[dPLP,d.y] 

|/lH-|Jll<|/j + |J| |Ji| + |Ji|<|/| + |J| 


Y dPLP<^ld.^dPLP + 2 dPLP<^l[dPL-’fd.f[ 

Il^l2 = I,J-i + J2 = J Il+l2 = i 

|Jll + |Jll<|/| + |J| Ji+J2 = J 


Y dPLP<^ld.ydPLP + 2 

Il+l2 = I,Ji + J2 = J Il+l2 = I 

|Jll + |Jll<|/| + |J| Ji+J 2 = J 


^ {dPPP<!>],) ePf,dsdPL-^i 


Hence, linear combinations of and {dPL'^'-^'p)9Pj, and Ji + J 2 = J, which 

yields (ILSI) . Note that 9'Pj, are constants, so that 

dPL-P{dPLP<^l9Pj,^) = 9-!^f^dPL'PdPLP^]i. 

By definition, <I>^ is a homogeneous function of degree zero, so that is again homogeneous 

but with degree ^ 0. We thus arrive at (HU). □ 


Lemma 4.4 (Algebraic decomposition of commutators. IV). Within the future cone 3C, the 
following identity holds 

(4.10) 2 ai^jd,L\ 

\J\<\i\ 

where the coefficients a(.j are smooth functions and satisfy (in %) 

(4.11) \d^^LPai]\ ^ C{\I\, I J|, |/i|, I Ji|)t-I^^l. 

Proof. This is also by induction. Again, when |/| = 1, (14.101) together with (14.111) are guaranteed 
by (|4.4I) . Assume that (14.101) and (I4.11|) hold for |/| ^ to, we will prove that they are valid for 
|/| = TO + 1. We take = LaL'^ with | J| = TO, and obtain 

[L^dJ =[KL\d,] = La{[L\d,]) + [La,d,]LJ 

=La( Y 

= 2 L,AJ,d,A'+ Y Aj'Lad,A' 

|,7'1<|J1 1.7'|<1.71 
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SO that 


[i^^c]= X 


+ X 


X 






t,7'l<|Jl 


- X 


+ X 


X 

Aj'Q-ibldL^' +QldbL^- 



\-T\<\A 


1,7'|<1.7| 



In each term the coefficients are homogeneous of degree 0 (by applying (j4.11|) b and the desired 
result is proven. □ 

The following result is also checked by induction along the same lines as above, and so its proof 
is omitted. 

Lemma 4.5 (Algebraic decomposition of commutators. V). Within the future cone 3C, the fol¬ 
lowing identity holds: 

(4.12) [d^^dJ=t-^Yi 

where the coefficients p^j are smooth functions and satisfy (in 3C) 

(4.13) \d^-L^-pij\ ^ C{\I\fJUhUJf)t-\^^\. 

4.2. Estimates for the commutators. The following statements are now immediate in view of 
(14.51) . (14.6[) . and (|4.10() . and (|4.12() . 

Proposition 4.6 (Estimates on commutators. I). For all sufficiently regular functions u defined 
in the future cone OC, the following estimates hold: 

(4.14) \[d^LJ,da,]u\^Ci\Il\J\) 1] \dpd^LJ'u\, 

\J’\<\J\,P 


(4.15) 


\[d^L\d,]u\^C{\I\,\J\)( 2 \d,d^'L^'u\ + t-^Y. \d^'L^'u\\ 

\ |J'|<|J|,a IUS|C| / 

If'isiu 


(4.16) \[d^L^,d^u]\^C{\I\,\J\)t-^ dpd^'L^'u +C{\Il\J\) Y 


D'KUI 


3,|P|SIU 

|J'|<|J| 


(4.17) 


\[d^L-fdo.dp\u\^C{\I\,\J\) Y L-^'u\, 


|/|S|PUJ'|<|/| 


(4.18) 

\[d^L'’,d^dp\u\ + 




Ci\I\,\J\)( Xj 1 ^ X u\ + t ^ X {d.yd^L'^un. 

V C.^T.I/'ISIJ-I c,7,|/'|<:|/| 7,|/'|s;PI / 


|j''l<|.'l U'NUI 

Further estimates will be also needed, as now stated. 


U'KUI 


Proposition 4.7 ([Estimates on commutators. II). For all sufficiently regular functions u defined 
in the future cone 3C, the following estimate holds (for all I,J,a) 

(4.19) \d^L'^{{s/t)dau)\^\{s/t)dad^L'’u\+C{\I\,\J\) X \{s/t)dpd^ L-^ u\. 

jJ'IsDI 

Recall that the proof of the above result (given in [22]) relies on the following technical obser¬ 
vation, concerning products of first-order linear operators with homogeneous coefficients of order 
0 or 1. 
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Lemma 4.8. For all multi-indices I , the function 

:= {t/s)d^LJ{s/t), 

defined in the closed cone % = {|x| ^ i — 1}, is smooth and all of its derivatives (of any order) are 
hounded in JC. Furthermore, it is homogeneous of degree rj with ly ^ 0. 

5. Initialization of the bootstrap argument 

From this section onwards, we begin the proof of Theorem 11.11 which is a rather involved 
bootstrap argument along the lines of the method presented in [^. We fix some integer iV > 8 
throughout, and first summarize our strategy, as follows. 

Let (u, v) be the local-in-time solution to the Cauchy problem associated with the system (11.1111 . 
From a standard local existence result (cf., for instance, [22l Section 11]), we can construct a local¬ 
time solution from the data given on the initial hypersurface and, consequently, guarantee that on 
the initial hyperboloid and for all |/| -I- | J| ^ TV, 

^ Coe, Em{2,d^L-^vy^^ ^ Cos 

for some uniform constant Co > 0. On some (hyperbolic) time interval [2, si], we can thus assume 
the following energy conditions for some constants Ci,£,6 > 0 (yet to be determined): 


(5.1) 


Em{s,d^L\y/^ ^ Cies'=^ 

\J\ = k, 

\I\ -1- 1 J| ^ TV, 

wave / high-order, 

Em{s,d^L-’uy^^ iiCie, 


|/| + |j|^TV-4, 

wave / low-order. 

Em{s,d^L'’vy^^ ^ Cies^/^+fci^ 

\J\ = k, 

\I\ + \J\ ^ TV, 

Klein-Gordon / high-order. 

Em{s,d^L-’vy^^ ^ Ci£S^^, 

\J\ = k, 

|/| + |J| < A-4 

Klein-Gordon / low-order. 

We will prove that on the same interval the following improved energy bounds are valid when e is 

sufficiently small and Ci > Co with ^ (5 ^ 

^ (fixed once for all): 


(5.2) 

Em{s,d^L-^uy/^ ss ^Cies'=^ 



\J\ = k, 

|J| + 1 J| ^ TV, 

wave / high-order, 



|/| + |J|^TV-4, 

wave / low-order. 


\J\ = k, 

|J| -f 1 J| < TV, 

Klein-Gordon / high-order. 


\J\ = k, 

|/| + \J\^N-4: 

Klein-Gordon / low-order. 


Once this property is proven, we set 

Si := sup |s /(EH) holds on [2, s]| 


and we can deduce that si = -l-oo. Indeed, by a continuity argument, Ci > Cq implies si > 2. 
Again by a continuity argument, we deduce that when s = si, at least one of the inequalities (EB 
must be an equality. But, when (15.21) holds, none of them can become an equality. This means 
that Si = - 1-00 and the rest of our work consists of proving (lOD . 

Proposition 5.1 (Formulation of the bootstrap argument). Given any integer TV > 8 and < 
i5 < there exist constants Ci,£ > 0 satisfying eCi < 1 such that any local-in-time solution 
{u,v) to (Hii), defined in the time interval [2, si] and satisfying the energy conditions (15.11) for 
some ee (0,eo]) o,lso satisfies the improved energy bounds (15.2p . 

The remaining text is devoted to the proof of this proposition, which we decompose into three 
parts. First, we derive a series of and sup-norm estimates directly from (EB, and from the 
Sobolev inequality on hyperboloids (12.141) and the commutator estimates (i.e. Pronositions ld.Gl and 
1471) . Second, we improve the sup-norm estimates by using (EB and EB- Finally, we combine 
the improved sup-norm estimates and the estimates established in the first part and we get 
the improved energy estimates ()5.2I) . 
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From now on, we consider a solution satisfying the energy bound (15.ip and throughout 
we set I J| = fc. 


6. Basic estimates 

6.1. Basic estimates of the first generation. Throughout, we always assume that 2 ^ s ^ 
si. We begin by stating the type estimates provided to us by the energy assumption (15.ip . For 
all |/| + I J| ^ with \J\ = k, we have the high-order bounds 

~ CiSS^^, 

(6.1) 

where the last estimate implies, for all |/| -f | J| ^ — 1 and | J| = k, the following estimate 

(6.2) 

as well as, for | J| + | J| ^ A^ — 4 with | J| = k, the low-order energy bounds imply: 

||(s/t)d„d^L“'u||i2(3^^) -t ||(s/^)d„d^L'^u||i2(j<;^) < Cie, 

(6.3) ||(s/t)dad'^L‘^u||i2(3^^) -I- ||(s/^)d„d‘^L“'w||i2(j^^) < C'les'"^, 

~ CiES^^, 

In addition, they also imply, for all |/| -I- | J| ^ A^ — 5 with | J| = k, 

(6.4) < Cis'=^ 

6.2. Basic estimates of the second generation. The following estimates are obtained 
by applying the above energy estimate combined with the commutator estimates presented in 
Proposition 14.61 For all |/| -I- | J| ^ A^ with | J| = k, we have the high-order bounds 

\\{s/t)d^L'^dau\\i^ 2 (^^^) + ||(s/^)d'^T‘^d„u||i 2 (j^^) < Cies^^, 

4- \\d^L''dj_u\\i^2^j^^f < Ciss'"^, 

(6.5) \\{s/t)d^L^dav\\L2(^^^-^ + ||(sA)d^T‘^d„u||i2(„^) < 

which, for |/| -f | J| ^ A^ — 1 with | J| = k, imply the low-order bounds (for instance, by expressing 
= La in the first inequality): 

\\d^L'^dav\\i^2(^j^^-f + \\d^L'^d^v\\i^2(^^^'f < , 

\\{s/t)d^L^dad(iv\\L 2 ^:}{^) + \\{s/t)d^L-^d^dpv\\L 2 ^j{^) < CiesL'^+’^^, 

\\sd^L’^d^di,v\\i^2(^j^^'f + \\sd^L’^ 


( 6 . 6 ) 
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Observe especially that, for derivation of the last term, we write 

= d^L-^{t~^Lad^v) = ^ L’’^ La^vd^^ L'’’^ {t~^) 

Il+l2=I 
Jl +J2 = J 

and then we proceed by homogeneity (by noting that is homogeneous of degree —1 and its 
derivatives are homogeneous of degree ^ — 1): 

\\sd^ 2 \\st-^d^^L^^Lad^v\\L2(^^) < ^ 

Uiisiii iHisin 

UiKUl OilsUI 


We also have, for |/| + | J| ^ — 2, 

(6.7) Wd^D’dML^.i^s) ^ C'i£s1/2+'=^. 

For |/| + I J| ^ — 4 with | J| = fc, we have 

||(s/t)d^L“'daM||i2(3^^) + \\{s/t)d^L-^d^u\\L2(^^) < Cic, 

L '’< Cie, 

(6.8) ||(s/t)d'^L'^d„i;||i 2 ( 3 ^^) + \\{,s/t)d^L'’< Ci£s^\ 

ll<5"^7,‘^d_i_u||^2(3^^) + < Cies^^, 

~ C'les^'*. 

For |/| + I J| ^ A^ — 5, I J| = fc, we have 

+ Wh^d^L^^K.) ^ Cies('=+i)^ 

||d^L^d„i;|U2^) + ||d^L^d„u|U2^) < Cics'=', 

(6.9) ^ ^ 

||(s/t)d^L'^dad/jw||i2(j^;^) + ||(s/t)d^L-^d„d^u||i2(j<;^) < C'ies'=^ 

\\sd ^+ \\sd^ . 
For |/| + I J| ^ N — 6, we have 


\\d^L'^dadfjv\\L‘j(D<^) + \\d^L'^d^dpv\\L2^:K,) ^ C'les'=^ 


6.3. Basic sup-norm estimates of the first generation. We combine the Sobolev inequality 
on hyperboloids (I2.14p with our estimates. In view of the high-order bounds, for |/| -I-1 J| ^ 
N — 2 with I J| = fc we have 


( 6 . 11 ) 


sup {t^^^s\dc.d^L\\) + sup {t^^^s\dJ^L\\) < Ci£s^’^+^'>\ 
sup L'^u\) + sup (t^^^|d_i_d^L‘^u|) < 

sup {t^^^s\dad^L-^v\) + sup {t^^^s\dJ^L^v\) < Cies^^^+^^+^'>\ 
sup + sup , 

J<s 0<s 


sup(t^/^|d^L'^z;|) < C'ies^/2+('=+2)^. 


For |/| 4- I J| ^ N — 3 with | J| = k, we have 

(6.12) snp(t^^^\do,d^L-^v\) < . 

Ms 
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(6.13) 


From the low-order bounds, for |/| -I- | J| ^ — 6 with | J| = k we have 

sup {t^^'^s\dad^L'^u\) + sup s\^d^L'^u\) < Cie, 

■K, -K, 

sup {t^/‘^\dJ^L^u\) + sup < Cie, 

Ms 

sup -t sup 

Ms Ms 

sup(A2|d_L<5YY|) -bsup(A2|d„(5Y'^u|) < 

Ms Ms 

sup(t^/^ld^L-^vl) < C'ies('=+2)‘5. 

Ms 

For |/| -I- I J| ^ — 7 with | J| = k, we have 

(6.14) sup (A^ldadY'^ul) < C'ies('=+2)‘^. 


Ms 


6.4. Basic sup-norm estimates of the second generation. For |/| -I- 1 J| ^ N — 2 with | J| 
we have the high-order bounds 

sup (t^^^sld^L-^dc^ul) + sup (t^^^sld^L-^d^ul) < C'ies('=+2)‘5, 

Ms Ms 

sup + sup {t^^‘^\d^L-^d_^u\) < , 

:Ks o<s 

(6.15) sup {f^‘^s\d^L'^dav\) + sup (t^/^s|<5YY„u|) < C'iesi/ 2 -Kfe-H 2 ) 5 ^ 

3~Cs 0~Cs 

sup(A2|dYY,u|) < Ciesi/2+(fe+2)5^ 

sup(A2|dY^u|) <Ciesl/2+(fe+2)6 

Ms 

and, for |/| -I- | J| ^ A^ — 3 with | J| = k, 

supft^/'Id.dY-'ul) < Ciesi/"+('=+3)^ 

Ms 

sup(t’5/2|(5YY„u|) < C'iesi/2+('=+3)^ 

Ms 

(6.16) sup -t sup (A2|dYY„u|) < C'iesi/2-Kfe-H2)5^ 

sup{t^^^s\d^L-^dad0v\)+snp{f^^s\d^L‘^d^dov\) < Cis^/2+(fe+2)5^ 
sup (A^sIdYY^dfcul) -t sup < C'lS^/^-Kfe-HS)^^ 

Ms Ms 

For |/| -I- I J| ^ A^ — 6 with | J| = fc, we have 

sup (t^^'^s\d^L'^dau\) + sup s\d^L'^d^u\) < Cie, 

Ms Ms 

sup L'^d^u\) + sup (t^^^|d‘^LY_i_'u|) < Cie, 

3<s Ms 

sup s\d^L'^daV^) + sup 

Ms Ms 

sup(A2|(5YY|) < C'ies('=+2)'5. 

<TJ- ^ 


(6.17) 
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In addition, for |/| + | J| ^ IV — 7 with | J| = k, we have 

sup{t^^^\d^L'^d^v\) < C'les('=+3)^ 
(6.18) sup + sup d^v\) < Ci£s('=+^)^ 

'' ' T-/" ‘TT 

DR q g 

sup (t^'^^s|d^L'^dQ,(5/3u|) + sup s\d^L'^d^dav\) < , 

sup {t^^^s\d^L-^d^di,v\) + sup s\d^d^d.v\) < Cis^'^+^'>\ 

TC, 

Moreover, |/| + | J| < iV — 8 with \ J\ = k, we have 


sup (t^^'^\d^L‘^dadpv\) + sup L'^< C'is^*"'-2)<5^ 

(6 19) 

sup (t^^'^\dc,dpd^L-^v\) + sup < Cis^’^+'^'>^. 

6.5. Estimates based on Hardy’s inequality on hyperboloids. We now substitute the basic 
estimates in Hardy’s inequality (12.171) and find 

(6.20a) ||s < Cqe + Ciss^^, \L\ ^ N, 

as well as the inequality (which will not be used in the following) 

(6.20b) ||s“^L'^M|| 2 , 2 ( 5 r„) ^ (Co + Ci)e + Cielns, |L|^A^ —4. 

6.6. Estimate based on integration along radial rays. By the first estimate in (16.131) and 

since (in the domain of interest), we obtain 

\drd^L'^u{t, cc)! < — r)“^/^, |/| + | J| ^ iV — 6. 

Then we integrate this inequality in space along the rays (t, Aa;)|osjAsjt-i for any x e §^: 

(6.21) |d'^L'^u(t,a;)| < Ciet“^(t — r)^/^ ^ Ci£t“^/^s, |/| + | J| ^ — 6. 


7. Refined sup-norm estimates 

7.1. Overview of the analysis in this section. We now proceed by using the structure of the 
nonlinear wave system under consideration, and relying now on sharp sup-norm estimates. In the 
following sections we are going to establish the following estimates: For |/| ^ — 4, we have 

(7.1a) sup t|M| < Cie, 


(7.1b) sup ((s/t) -I- sup ((s/t) v\) < Cie, 

0<s Ks 

and, more generally, for |/| + | J| ^ A^ — 4 with | J| = k, 

(7.2a) supt|L‘^u| < Cies^"^, 

(7.2b) sup ((s/t)-3/2+45^3/2|^^^7^J^|^ g^p ((s/t)-l/2+4^t3/2|^/^J^|^ < C'les'^^ 

Ms 

The property dzu is essentially a special case of n -. we will establish it first and it will next 
serve in the proof of (17.2F done by induction on k. 

The sup-norm estimate for the Klein-Gordon component (13.71) and the sup-norm estimate for 
the wave equation (13.21) will now be used. We proceed with the following calculation: 

(7.3) - n{d^L^u) = P°‘f^d^L\do,vdpv) + Rd^L\v‘^), 

(7.4) 


n{d^L\) + H°‘^udadpd^L\ + c^d^L\ = [i7“'^udad/3, 
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We also recall by (1X71) . the Ri terms in this context (with h°‘^ = read as follows: 


Ri[d^L-^v] 

R2[d^R^v] 

Rsid^L-^v] 





5a5b + 


4sl/2 


^ Zj gl/2 


d^L\ 


- s^l‘^{2li'%'dbd^L-^v + T%dbd^L-^v + R-’v), 

^2a:“s^/^doda + 


Hence, the following four terms must be controlled: 

(7.5) L'' [davdpv), d^L''{v‘^), Ri[d^L^v^, [H°'^udadp,d^ L'^^v. 

7.2. First improvement of the sup-norm of the wave component. We now present esti¬ 
mates which use only the basic sup-norm estimates already established in Sections 16.31 and 16.41 
We hrst estimate the terms L'^ {davdpv) and d^L‘^{v^). 


Lemma 7.1. If the energy hounds dSU) hold, then for all |/| -I- | J| ^ N — 7 with | J| = k, the 
following estimate holds in the region 3C[2,si] ■’ 

(7.6) \d^L-^{do.vdpv)\ + \d ^^ C{Cieft-^s^^+‘^^^. 


Proof. We have 

(7.7) d^L-^{do,vdisv) = ^ R-^^v, 

Il+l2 = I 
Jl +J2 = J 

where, in the right-hand side, each term satisfies |Ji| + I/ 2 I = |/| and | Ji| + | J 2 I = | J|. Then we 
obtain 

where we have used the third inequality in (16.181) for each term. The estimate of R"^ (y^) is 
derived similarly. □ 


We improve the bound on u, as follows. 


Proposition 7.2 (First improvement of the sup-norm of the wave component). For |J| -I- | J| ^ 
N — 7 one has 


(7.8) |d^L'^u(t,a:)| < -f {Cief{s/t)^'^+^^h-h^'^+^'>^. 


Proof. The proof is a direct application of the sup-norm estimate for the wave equation (13.211 . 
First of all, R'’u solves the Cauchy problem 

nd^R\ = d^R''{P°‘f^davdpv) + d^R'’{Rv^), 

R‘^u(2,x) = Uo{I,J,x), dtd^ R'^u(2,x) = Rfi{I,J,x), 

where Uo{I,J,x) and Ui{I,J,x) are restrictions of d^R'^u and d^R'^u on the initial hyperplane 
{t = 2}. We remark that they are linear combinations of d^,, u and dtdl. u with |/'| ^ |/| -I- | J|. 
Hence, u is decomposed as follows: 


u(t, x) = Wi{t, x) + W 2 {t, x) 

with 

□L'^zi ;2 = R\P^^do.vdpv) + R^{Rv'^), 
W2(2,x) = dtW2(2,x) = 0, 

while 

□wi = 0, 

wi(2,a;) = Ui{I,J,x), dtWi{ 2 ,x) = 1/2(1, J,x). 
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The sup-norm bound for wi comes directly from the explicit expression of the solutions (cf., 
for instance, [32]) while for W 2 we apply (jUj). Observe that for the terms in the right-hand side: 

\d^L-^{P°‘^dc,vdfiv)\ + \d^L^{Rv‘^) \ 

<(C'^£)2^-2-(l-(2 + fe/2)5) _ ^^-l + (l + (2+fe/2)5) 

in 3C[2 ,s] c 3<i[2_si]- Recall that this estimate also holds in {{t,x)\r < t — 1,^ sf,t > 2}. 
Then by (13.21) . the desired result is proven. □ 

We next estimate the terms Ri[d^L'^v]. 

Lemma 7.3. For |/| -)- |J| ^ fV — 4 with |J| = k, the following estimates hold in 3C[2,si]' 

(7.10a) \Ri[d^L-^v]\ < Cie(s/t)3/2g-3/2+(fe+4)5^ 

(7.10b) \R2[d^L-^v]\ < C'ie|tu|(s/t)3/2g-3/2+(fe+3)5 (^^^)2(^/^^3/2g-3/2+(fe+3)5^ 

(7.10c) \R 3 [d^L-^v] \ < {Cief{s/t)s-^+^'^+^^^ + Cis\tu\{s/tf^^s-^^^+^'^+^'>^. 

Proof. The proof is a substitution of the basic sup-norm estimates into the corresponding expres¬ 
sion. We begin with Ri and focus first on dadbd^L'’v: 

dadbd^L'^v =t~^La{t~^Lbd^ L'^v) = t~^ La{t~^d^ LbL'^v + t~^[Lb,d^]L'^v) 

= t~^La{t~^)d^LbL'^v + t~‘^Lad^LbL’^v + t~^La{t~^) [Lb, d^]L‘^v 

(7.11) +t-^La[Lb,d^]L-^v 

= t-^La{t-^)d^LbL\ + t-^d^LaLbL-^v + t-'^[La,d^]LbL\ 

+ t~^La{t~^)[Lb,d^]L-’v + t~‘^La[Lb,d^]L'’v. 

For the last term, we apply (14.51) as follows: 

t~‘^La[Lb,d^]L-’v =-t~‘^ ^ \lj,Lad^'L-^v 

= -t-2 2 \ij,d^'LaL\-t-^ Y, Ki'[La,d^']LJv 

b'i=shi b'i=sbi 

= -t"2 Y \ird^'LaL-^v + r^ ^ XIp ^ L-^v, 

\i'\m h"i«h'i 

and the term [La,d^]LbL'’v is bounded in the same manner. Then we conclude that 

\t-^La[Lb,d^]L-^v\iiCt-‘^ 2 \d^'L-^'v\. 

|j'is|j|+i 

In view of the inequality |Lo(t“^)| ^ Ct~^ (in 3C), the terms in the right-hand side of (17.111) are 
bounded by Ct~^Yi i^'isi^i ^|- Then, by the last equation in (16.111) . we have 

|J'ISDI +2 

\s^^^dadbd^L^v\ < C'iC(sA)7/2^-3/2+('=+4)^ 

and, similarly, 

\x^s-^^^dad^L-^v\ < C'le(s/^)3/2s"^/^+('=+3)^ 

\x‘^x'’s-^^^dadbd^L\\ < C'ie(s/t)3/2s-3/2+('=+4)^ 

\s-^/^d^LJv\ < Cie(s/t)3/2s-3/2+('=+2)5, 

So we conclude that 

|Ri[d'^L“'u]| < < Cie{s/tf^'^s-^^'^+^'‘+'^'>^. 

For the derivation in the paragraph above, let us provide some more details by observing that, 
for the first term, 

[dadbd^L'^v] = \t~^La{t~^Lbd^L-^v\ < t~‘^\LaLbd^L'^v\ + t“^|La(t“^)| • \Lbd^L'’v\. 
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Here, we remark that 

S^/h-^\LaLt,d^L^v\ < t- 2 Cies 3 / 2 i- 2 gl /2 + (|/|+2+2)5^-3/2 < ^^^^^7/2^-3/2+(/c+4)5_ 

In the first inequality, we have used the fact that 

\LaLbd^L\\ < \d^LaLbL\\ + \[LaLb,d^]L\\. 

For the first term in the right-hand-side of the above inequality, we get the upper bound 

(for I J| -t I J| -f 2 ^ — 2) and for the second we recall the estimate on commutator (using (14.51) 

twice), and see that 

\[LaLb,d^]L-^v\ < ^ d^'\L-^'d^'L-^'v\ < . 

|J'K|J+i| 


The estimates of i ?2 and Rz are quite similar. We just need to observe that, by (j6.21l) and by 


ittOO 


.-^aO 


-^ab I 


recalling that \H \ ^ C{t/s)‘^, \H \ ^ C{t/s), and \H | ^ C, we obtain 


(7.12) 

and 


\h I = |i7 u| ^ |/i° I < Cist \h°' u\ < Ciet 

Idodad^L-^v] ={s/t)\dtdad^L-^v\ 

^{s/t)t~‘^\Lad^ L'^v\ + t~^{s/t)\dtLad^ L'^v\. 

As was done in (17.lip and by applying (14.51) and the fifth equation in (16.151) we find 
(7.13) \dodad^L-^v\ < 

Equipped with (17.121) and (I7.13p . we see that in R 2 [d^L'^v], 

< (C'ie)2(sA)s-"+('=+2)^ < (C'ie)2(sA)3/2s-3/2+('=+2)^, 

< C'ie|tu|(s/t)3/2s-3/2+('=+2)'5, 

< (C'ie)2(s/t)5s-2+('=+4)5^ 

< C'ie|tu|(s/t)3/2s-3/2+('=+2)'5, 

= 0 , 

while, in the expression R 3 [d^L‘^]v, 

\Tf\^s^^^dodad^L-^v\ < C'le|^^t|(s/^)3/2s-3/2+('=+3)^ 

< (Ci£)"(s/t)3/2s-3/2 + (fc+3)5^ 


\s^^‘^h°^dod^ L'^v 
s^^^lTf’dodbd^ L-^v 
s^^^\h‘^^dadbd^ L'^v 


dod^L\ 


^3/20 IjJy 


s-^/^\h^\‘^x^dadbd^L-^v\ < (C'ie) 2 (s/t)s- 2 +('=+ 4 )< 5 . 


□ 


7.3. Second improvement on the wave component and first improvement on the Klein- 
Gordon component. We now establish ()7.1ap - (l7.1bp and, for latter use, we first derive the 
following improved estimates on the terms Ri. 

Lemma 7.4. For |/| -I- |J| ^ A^ — 4, | J| = k, the following estimates hold in lK[ 2 ,si]- 
(7.14) ^ \R,[d^L-^v]\ <C'ie(s/f)3/2s-3/2+('=+^)^ 

i=l 

Proof. This is a combination of Lemma 17.31 and (17.81) (take k = 0 then considering the condition 
Ci£ ^ 1) and the fact that in X, ^ s ^ f. □ 

Then we need to estimate the term |/i( a;(A)| in Proposition 13.31 
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Lemma 7.5. The following estimate holds for (t^x) e 3(!|-2,si]' 


(7.15) 


f \KJx)\dx 

Sn 


<Cie, 


1-00 


-rrOO , 


where ht^x{X) := h {Xt/s, Xx/s) = H {Xt/s, Xx/s)u{Xt/s, Xx/s). 

Proof. With the notation of Proposition 13.31 we have = H^^u, and we observe that 




-ttOO, 


tOO, 


tOO , 


Note that H {Xt/s, Xx/s) = H {t, x), so that H is constant along the segment (At/s, Xx/s), sq ^ 
A ^ s. So we find 

K.xiX) = H°°{t,x){t/s)dj_u{Xt/s,Xx/s), 

and we conclude that 

l^t.x(^)l < C{t/sf\dj_u{Xt/s,Xx/s)\. 

Next, we observe the identity 


(7.16) 


S^ x^ s^ x^ 

dj_u = -^StU + —^aU = -^3tU -f -^LaU 


and, by the first inequality in (16.171) and (17.81) with L'’ = La, 

|d_|_M| < Cie{s/t)t~^/'^ + Ci£{s/if‘^t~‘^s^^. 

Therefore, we obtain 

\K^,,{X)\ < Ci£{s/t)-^^'^X-^^^ + Ci£{s/t)-^+^^X-^+^^. 

Then, to apply the sup-norm estimate for the Klein-Gordon equation (1^ . we proceed as 
follows. In the range 0 ^ r/t ^ 3/5, we have 4/5 ^ s/t ^ 1, and 

r |/i/^(A)|dA < Cie r A-3/2 dX < Cis. 

JsQ ’ J2 

In the range 3/5 < r/t < 1, we obtain 

r |/r/^(A)|dA <C'ie(s/t)-i/2 r A'^/^ dA + C'ie(s/t)-i+'^'^ f X-‘^+^^ dX 

J Sq Jsq *JSo 

< C'ie(s/f)“^/^S(/^^^ -I- C'ie(s/t)“^+^'^SQ 
— . / t+r > 


We recall that, when 3/5 < r/t < I, sq = ... 
desired result is established. 


t/s, so that \h/a;{X)\dX < Cie, and the 

□ 


Now we give a second application of the sup-norm estimate for the Klein-Gordon component 

(ITTI) . 

Proposition 7.6 (Second improvement on the wave component and first improvement on the 
Klein-Gordon component). The following estimate also holds in lJC[ 2 ,si]' 

(7.17a) |u(t,a;)|-l—|d_|_u(t,x)| < Cie(s/t)^“^‘^s“^/^, 

s 

(7.17b) |M(t,x)| < Ci£t~^. 

Proof. We rely here on (13.71) and the sup-norm estimate for the wave equation (13.21) , and we first 
establish (I7.17al) . In view of (17.151) . we have 

Jsq ^Sq 

\K,xis)\e^^^‘'ds < Ci£. 

'Jsq 
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On the other hand, to estimate F{s), we write 

F{s) = r (i?i[n] + i? 2 [n] + -R3[i^])(At/s, Aa;/s) dA 
Jso 

< Cie (s/t)3/2A-3/2+7^dA < 

'JSq 

Now for 0 ^ r/t ^ 3/5 we see that 4/5 ^ s/t ^ 1 and sq = 2, and we have 

F{s) < Cie(s/t)3/^So < Cie{s/tf-'^^. 

For 3/5 < r/t < 1, we see that sq = a/ < t/s, so that 


(7.18) 


F(s) < Cie(s/t) 


2-7(5 


Then, by combining (I7.15|) . (17.181) and (13.7L we conclude that (I7.17al) holds. On the other hand, 
(I7.17bp follows directly from substituting (|7.17all into (13.21) . 

Let us explain in more detail the above argument. In the equation 

□u = PafidaVdpV + Rv^, 

we need to estimate \P°‘^davdpv\ and First we rewrite the expresison P^^^davdpv in the 

semi-hyperboloidal frame: 

P'^^da.vdpv =P°‘'^d^vdpv 

=^^dtvdtv + P/'^d^vdtv + P^^dtv^v + P/^^^vd^,v. 

The last three terms in the right-hand-side can be controlled by applying the first and the third 
inequalities in (I6.18P 

\P’^°^d^vdpv\ + ^ C'(C'ie)^t"^s®^. 

For the first term, we see that 


s‘ 


dtV = — [ dj_v • 


-laV 


Then by (I7.17al) and the third inequality in (I6.18p . we obtain 

\dtv\ 

This leads to 

\P°^dtvdtv\ ^ C{Cief{t - r)-^Hl/2-7S/2)^-2-il/2-7S/2)_ 
The term is bounded directly by (l7.17aF and we have 




+-3 


Then by applying (IQ) . the desired bound on u is guaranteed. 


□ 


With (iTlTbl) . we can improve again the estimate on Ri. Namely, the proof of the following 
estimate is immediate by substituting (I7.17bl) into (17.7|) . and using (13.2F 

Lemma 7.7. The following estimates hold: 

(7.19) 2 Ri[d^L\] < C'ie(s/t)3/2s"3/2+('=+'‘^^. 

2=1 
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7.4. Second improvement on the Klein-Gordon component. We now establish (I7.1bl) and, 
to do so, our first task is to estimate the commutator First of all, from the 

following identities 

y-2 j.™ci 

(7.20) d, = 5a = -^5^ + + d„ 

the following estimates are immediate: 

^ (t/s)^|d_i_d^L'^u| + (t/s)^^ \d_^d^L'^v\, 

(7.21) 

a 

Based on this result, we estimate the commutator [H°‘^udadi 3 ,d^]. (In the statement below, as 
usual, a sum over the empty set vanishes.) 


Lemma 7.8. The following estimates are valid in X for |/| + | J| ^ — 4 and | J| = k: 


(7.22) 


|[iJ“^uda5/3, d^L'^]v\ 

< Ciet-\s/t)-'^ 2 ^ \L-^^u\ \d^L'’^dc.di3v\ 

|/2l^UI./3 Jl + J2 = J 

+ Ciet-3/2(s/t)-3 ^ \d_^di3d^^L-^^v\ 

|/2l + |J2l^UI + |J|-l,/3 
N-7^\l2\ + \J2\^^-^ 

+ {Cief{s/tf^^s-^+^’‘+^'>^. 


Proof. We write the decomposition 

[H‘^>^ud^dp,d^L'^]v= ^ H^^d^^d-’^ud^^L-^^dadpv + H°‘'^ud\[L^,dadp]v) 


Il+l2=I 
+J2 = J 
irp + Dilsi 


= :To + T7. 

We recall that by (14.171) and (I7.17bp . Tr is bounded as follows: 


\T7\<Ciet-^ Y, \d^da,dpL 




a:,/3 

1^21^1^1-1 


<Ciet ^{t/sY ^ |d_i_d/3d^L'^^u| + Ciet ^it/sY ^ \d_^dpd^L'^'^v\. 

P a,13 

J 2 ISDI -1 UalsDI-i 

The second term in the right-hand side is bounded as follows: 

Cist-\t/sf\d^di3d^L-^^v\ < C'ier^(t/s)2C'ier®/2si/2+(fe+3)5 

<(C'ie)2(s/t)3/2s-3+('=+3)<5. 

We then write 

\To\^ Y \d^^LJ^u\\d^-LJ-da.dpv\+ Y \L^^u\\d^LJ-d^dpv\ 


Ji + /2 = -7,|Ji|>1 
Ji + J2 = J,c.,/3,7 


Jl + J2 = 'J 


= :Ti + T 2 . 

Then we see that Ti is again decomposed as follows: 

Ti= ^ \d^^L''^u\ Id^^L-^^dc^dpv] + ^ \d^^L-^^u\ \d^^L'^^dad0v\ 


/l+/2 = J,|/llSl 
Ji + J2 = Ac«,/3,7 

T2l + D2ls:«-8 


Ji+J2 = J,c<,/3,7 
N-7<S|/2l + D2l«'^-5 


=:T3 -f n. 


We have 
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where we applied (16.191) and the fourth estimate in (16.181) . Then by applying (17.211) and in view 
of (|4.17|) . the term is bounded by 

r 4 ^ Y, \d^^L-^^u\\dadpd^^L-’^v\ 

|4ll + |42l«l4|,c«,/3,7 
^V-7s:|/2l + l J 2 NW -5 

Xj \dj_d/3d^'^L'^^v\ 

|4ll + |42lS|J|<»,/3,7 
N-75:|/2l + |42|s:iV-5 

+ {t/s)^ X \d^^L-^^u\\^d^d^^L^^v\=: Ts + Tg. 

/l+J2 = J‘,PllSl 
|4ll + |J'2lS|4|,<i,/3,7 

N-75:|/2l + |42lSN-5 

Then, in the expression T 5 , N — 7 ^ 1^21 + |•^ 2 | ^ N —5 implies |/i| + |Ji| ^ 3 ^ TV — 6 and recall 
|/i| ^ 1, so we see \d^^L'^^u\ ^ were \I[\ + |Ji| ^ 2 ^ — 6 . Then by the first 

estimate in (16.131) . we find 

T 5 < Ciet-3/2(s/t)-3 ^ \d^dpd^^L-^^v\. 

|72l + |J2l«|J| + |J|-l 
N-7SI/2I + I42ISW-5 

Furthermore, we have 

To < < (Cie)2(sA)s-"/2+(fe+2)5 

<(Cie)2(s/t)3/2g-3 + (fc+3)5 

and the desired estimate is established. □ 


We are now in a position to establish the desired bound (I7.1bl) . 

Proposition 7.9 (Second improvement on the Klein-Gordon component). The following estimate 
holds in 3C[2_sj] fof |/| < TV — 4: 

(7.23a) \d_id^v{t,x)\ < Cie{s/t)^^‘^~^^t~^^^ , 

(7.23b) \d^v{t,x)\ < Cieis/tf/'^-^H-^l^. 

Proof. We first discuss the case where |/| — 1 > — 7 and, in this case, using (17.221) 

\[H°‘^ud^di3,d^]v\ <Ciet-^^\s/t)-^ X \d_Ld/sd^^v\ + {Cief{s/tf^‘^s-^+^^ . 

\l2\^\I\-l.f3 

N-7s;\l2\s;N-5 

For all s G [soj s], using (17.191) and the above estimate we have 

F{s) ^ X f \Ri[^^v]{Xt/s, Xx/s)\dX + r X^^^\[H°‘^udcxdi 3 ,d^]v\dX 

i=l'lso 9so 

<Cie(s/t)3/2 r dX + {Cxef{s/tfl‘^ { X-^/'^+'^^dX 

^Sq Jsq 

+ Cie{s/f)~^^‘^ X f \§.x5pd^'^v{Xt/s,Xx/s)\dX 

172lS|7|-l./3 -^*0 


N-7^\l2\^N-Z 


<Ci£{s/t) 3/2 ^ r |d_i_d,ad‘^^ 7 ;(At/s, Ax/s)|dA + C'ie(; 

/oISl/l-lA -^*0 


ft) 


3/2 - 1 / 2+45 


/2ISUI-1.3 
Ar_7^| j2|^Al-5 


Case I: 3/5 < r/t < 1. In this case, sq = \ jif- ^ t/s and we have 


F{s) <Ci£{s/t)‘^-'^^+ Ci£{s/t)-^/'^ X f 1 ^ 15 / 3 < 5 ^"?^(At/s, Ax/s)|dA. 


72lS|l|-1.3 

iV-7s:|l2|5:N-5 
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We define 

and find 
(7.24) 


Vt.x(A) := ^ \d^di3d^'^v{\t/s,\x/s)\ 


N-7^\l2\^N-^ 


F{s) < Cie{s/tf-^^ + Cie [ dX, sq 

Jsn 


^ S ^ S. 


Recalling the sup-norm estimate for the Klein-Gordon component (1X71) in the case 1 > r/t > 
3/5, we obtain 

^ Jso ’ ^ 

We replace {t,x) by (At/s, Ax/s) with sq < A ^ s, we see that (At/s, Ax/s) is again contained in 
^[ 2 ,si] ■ Then (13.71) still holds, and so 

|d_Ld^x(At/s, Ax/s)| ^ C'(s/t)A"^/^(^F(A) + f J^(s)|ti; ,,(s)|eS»- 

J Sq 

which implies 

nX 

{Xt/sf^^\d_^d^v{Xt/s,Xx/s)\ < C'ie(s/t)"^/GF(A) + J^(s)|ti/^(s)|eS^ ds). 

^ Jso ^ 

Recall that (I7.15|) holds for 1 > r/t > 3/5 and that F is increasing, then 

r J^(s)|ti/^(s)|e5^ ^ F(A) f |h/,,(s)|eS»-ds < C'ieF(A). 

'JSq Jsq 

So we see that 

(At/s)^/^|d_|_d^r(At/s, Ax/s)| < Ci£(s/t)“^/^F(A) 
and, in combination with (j7.24l) , 

(At/s)^/^|d_|_d^x(At/s, Ax/s)| < C'ie(s/t)^/^“^'^-I-Cie(s/t)“^/^ f s~^^‘^Vt,x{s)ds, 

'JSQ 

which implies (by taking sum over N — 6 ^ |/| — 4): 

(7.25) Vt,x{X) < Cic(s/t)3/2-45 + c'ie(s/t)-i/2 f s-^/^Vt,x{s)ds. 

'JSQ 

Then, by Gronwall lemma, we see that 

s-3/Vi.,(s)ds < Cie(s/t)3/2-4^ p --3/2^CCMs/t)--'^ s; 9-^'^d9 

'JSQ 'Jsq 

< C'ie(s/t)3/2-^‘^ f 

Jsc 


--3/2 


< 




Here we recall that sq = y fry ^ i/s, then 

(7.26) yt,,(A) < Cie(s/t)3/2-4^. 

Now we substitute (17.261) into (17.241) . and obtain 

(7.27) F(s) < Ci£(s/t)^"^^ So ^ s ^ s. 

Then we apply the sup-norm estimate (13.71) in the case 1 > r/t > 3/5 and considering (I7.15p . 

(7.28) a^)| ^ Gi£(s/t)^“"^'^s“^/^, |d‘^r(t,x)| < C'ie(s/t)^“'^^s“^/^. 
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Case II: 0 ^ r/t ^ 3/5. In this case, 4/5 ^ s/t ^ 1 and Sq = 2, so the discussion is simpler. We 
just remark that as in the former case, 

F(s) <C'ie(s/t)3/2s-i/2+6-5 + ^ f \s^d^sh^i^xt/s,Xx/s)\dX 


N-7^\l2\^N-5 


<Ci£ + CieJ A-3/Vt.,(A)dA. 

Then by the sup-norm estimate (13.71) (with 0 ^ r/t ^ 3/5), 

|d_Ld^p(t,x)| J 

+ t-3/2(^F(s) -f J 

Then similar to the former case, we get (recall (I7.15p l 

nX nX 

|(At/s)^/^d_i_d^p(t,x)| < (Co-f C'i)eCiE J s~^/‘^Vt^x{s)ds < Cie + C^e^ s~^/‘^Vt,x{s)ds, 

provided by Ci > Cq, which implies 

n\ 

I^,x(A) ^ CiE -f ^ ^^'^^t,x{s)ds 

Then, Gronwall lemma implies Vt_x(A) < Cie and, therefore, 


|d_Ld^p(t,x)| < Cirf-3/2 < . 


And again, as in the former case, we see that \d^v{t^x)\ < Ci(s/t)^/^ 

When I/| — 1 < TV — 7, we see that in this case 

\[H°‘^udadp,d^^v\ < (Cie)^(s/t)^/^s“^"'"'^'^. 

A direct application of the sup-norm estimate (Irfl) combined with (I7.19P will give the estimate 
on d^v and dad^v. Finally, combining these two cases, we see that the desired estimates are 
established. □ 


7.5. Third improvement on the wave and Klein-Gordon components. We now establish 
(EH), by combining the sup-norm estimate for the Klein-Gordon equation (1X71) and the sup-norm 
estimate for the wave equation (13.2L together with an additional bootstrap argument. 

Proposition 7.10 (Third improvement on the wave and Klein-Gordon components). There exist 
constants C, £2 > 0 (depending only on N ^ 8 and the structure of the model system (jl.lll) ) such 
that if the bootstrap assumption ED holds for £ ^ £2 and CiS ^ 1, then the following estimates 
also hold for all s e [2, si] and |/| -f | J| ^ TV — 4, | J| = k: 

(7.29a) sup (t|L'^u|) < Ci£s*^, 

(7.29b) sup {{s/t)-^+'^^s^^‘^\dj^d^L^v\) + sup {{s/t)-‘^+'^^s^^‘^\d^L-^v\) < Ci£s'=^ 

(7.29c) sup {{s/t)-^+'^^s^^‘^\dad^L‘^v\) < Ci£s'=^ 

Furthermore, we see that by the commutator estimates in Proposition ESI the following refined 
decay estimates are a direct consequence of (I7.29cl) : 

(7.30) sup ((s/T)^/^+^'^t3/^|d^L'^daP|) < Ci£s'=^ |/| -f | J| ^ A^ - 4, | J| = k, 

(7.31) sup ((s/T)^/2-H75^3/2|^/^J^^^^^A < CiEs’^^, \I\ -f | J| ^ A^ - 5,1 J| = fc. 
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(7.32) sup L^dcv\) < Cies^^, \I\ + \J\ ^ N - 5, \J\ = k, 


(7.33) sup ((s/t)-3/2+75^3/2|^/^j^^^^^Q < C'les'^^ |/| + | J| < iV - 6, \J\ = k. 


Proof. We proceed by induction on | J| and introduce the notation 

Vfc,o(A) := sup sup ((s/t)“^+^‘^s^/^|d^L'^u|), 

2s:ss:X.|J|s:fc 
|/| + 0|5:JV-4 

Vfc,i(A) := sup sup ((s/t)“^+’^‘^s^/^|dxd^i'^u|), 

2s:4S:X.|J'|s:fc 
\I\ + \J\^N ~4 

and, with |J| ^ k, Uk{X) := sup 2 s:ss:A supj£ (t\L‘^u\). To begin with, we observe that by (I7.23P 

|J|^fc _____ 

and (I7.17bl) . there exists a positive constant C determined by the structure of the system (II.lip 
such that 

Vb,o(A) < CiE, Vb,i(A) < CiE, C7o(A) < Cic, 

That is, (17.291) is proved in the case where fc = 0. 

Then we suppose that for allO^j l^fV — 5, there exists a (sufficient large) constant 

Ck-i depending only on the structure of the system (11.111) and a positive constant e'^_^ such that 
for all e ^ s'k-n 

(7.34) y,-o(s) < Ck-iC,es^', 1/,.i(s) ^ Ck-iC,es^^, U, ^ Ck-iC,es^' 

hold on [2, si] with Ck-i depending only on k and the structure of (jl.lip . Then we will prove 
that there exists a pair of positive constant (Cfc, e).) depending only on N and the structure of the 
model system (11.111) such that if (15.11) holds with e ^ e). and CiS ^ 1, then 

(7.35) Vk,ois) ^ CkCies>^\ VkAs) ^ CfCies'^^, Uk^CkCies^^. 


We rely on a bootstrap argument. First, we observe that on the initial hyperboloid IK 2 , there 
exists a positive constant Co > 0 such that 


(7.36) 


max sup((2/t) ^ Co feCiC, 

|/| + D|SAr_4 ^ ' 

IJISfc ^ 

max sup ((2/t)“^/^'''^'^f^/^|dj_d^L‘^i;|) ^ Co feCiC, 

|/| + DKAr_4 nrf ^ 


max sup (t|L‘^u|) ^ Co,feCie. 


We choose Ck > Co,fc and set 


32,k ■= supjs e [2,si] Vfe,o(s) ^ CfcCics'"'^, 

VkAs) ^ CkCies'^^ Uk ^ CkCies’^^y 

By continuity, we have S 2 ,fe > 2. We will prove that for all sufficiently large constant Ck > 
max{Co,fe, Cfe_i, 1} the following bounds hold on [2,S2,fc]: 

(7.37) Ffe,i(s) ^ iCfeCles'=^ Cfc ^ ic^Cies''^ 

for sufficiently small e. Once this is proven, we conclude that S 2 ,k = Si- Namely, proceeding by 
contradiction, we see that in the opposite case at S 2 ,k < si, at least one of the following conditions 
must hold: 

VkAs) = CkCies'^^, VkAs) = CkCies’^^, Uk = CkCies'^^ , 
which contradicts the improved estimates (17.371) . 

It remains to establish (j7.37l) and we derive first the following estimate for | J| + | J| ^ — 4, 

I J| = j ^ fc (again provided 2 ^ s ^ 32,k) 

(7.38a) \dAA\ < CfcCie(s/t)2-^'^s-3/2+^'^ 




















44 


PHILIPPE G. LeFLOCH AND YUE MA 


(7.38b) \dJ^L^v\ < CkCie{s/tf-'^^ 

(7.38c) \L-^u\ < CkCierh^^. 

The derivation of (I7.38al) is direct from the decay assumption (17.341) and the induction assumption 
(17.341) . while (I7.38bl) follows directly from (I7.21|) . the decay assumption (17.351) or the induction 
assumption (17.341) : 

\dad^L'^v\ < {s/t)~^\dj_d^L'^v\ + (s/t)“^^ \d^d^L‘^v\ 

a 

<C'feC'ie(s/t)-2(s/t)3-75g-3/2+j5 ^^^^^^^^-2^-5/2gl/2+(i+3)5 

<C'feCie(s/t)^-^'^S-3/2+7'^ + C'i£(s/t)^/2g-2+0-+3)5 

<CfeCie(sA)i-^^s-3/2+J'^, 

where the first equation in (16.161) was used. On the other hand, (I7.38cl) is also direct from (I7.34|) 
and ()7.35p . 

Then we need the following two estimates for |J| + | J| ^ — 4, | J| = k: 

(7.39) \d^L^{da,vd0v)\+ \d^L^{v^)\ < (CfcCie)2t-2-(l/2-75+fc5/2)(^_^)-l + (l/2-75+fcV2)^ 

(7.40) \[H°‘^udadp,d^L'’]v\ <{CkCief{s/tf-'^^. 

The estimate (|7.39l) follows directly from (17.381) . We see that 

Il+l 2 =I 
J-]^ J2 = J 

<C7feCi£(s/t)2-7'5s-3/2+l‘^il-5c'^C'^£(5/t)2-75^-3/2 + |J2|5 
^{CkCief{s - r)-l + (l/2-75+fe<5/2)^-2-(l/2-75-fc5/2) 

and 

\d^{^davdi3v)\ ^ ^ \d^'^L'^^davd^^L'^^di^v\ < (CfcCie)^(s/t)^“^'*^s“^"'"^^ 

Il+l 2 =l 
+ •^2 ~ 

,^(C7^C7l£)2^-2-(l/2-75-fe5/2) _ ^)-l + (l/2-75+fe5/2) ^ 

The estimate of (I7.40p is also direct by substituting ()7.38|) . We recall (|7.22() an write 
\[H°‘^uda,di3,d^L‘^]v\ 

< Cie(s/t)-2t-lCfce(s/t)3-7^S-3/2+fe^ + (CkS)^ ^ ^-1^|Jx|^(^/^)1-75^-3/2 + |J.|5 

l■7l| + |J2|s:|J| 

+ C7i£r3/2(s/t)-Yfe£(s/t)3-^'^S-3/2+fe5 (C7i£)2(s/t)3/2s-3+(fc+4)5 

< (C'feCi£)2(s/t)2-7'5s-5/2+fe5^ 

where we have assumed that Ck ^ Ci. 

Now we substitute (17.391) into (13.21) and find that (similar to the proof of Proposition [T^]) 

(7.41) < Co,fcCi£t-3/2 + {CkCief{s/t)'^^t-^s'^^, 
which is 


(7.42) Uk{s) < Co.kCie + {CkCiefs>^^. 

On the other hand, the estimate on and \dad^L'^v\ is a bit more difficult. We see that 

ps 3 ps 

F{s) ^ Ri[d^L'^v]{Xt/s, Xx/s) dX + A^/^|[i7“^MdQ,d/3, (5^L'^]u|(At/s, Ax/s) dA. 

J sq sn 
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By the sup-norm estimate (I3.7L in the region X n {3/5 < r/t < 1}, recall that sq > Ct/s, we 
can calculate each term in the right-hand side of the above inequality with that aid of (17.401) and 
(17.141) and find that 

\F{s)\ < 

Then, we apply the sup-norm estimate (Irfl) with (17.151) and by the same procedure in the proof 
of Proposition (17.61) . we conclude that when 3/5 < r/t < 1, 

(7.43a) \d_^d^L^v{t,x)\ < {Cie + {CkCief){s/tf-'^^, 

(7.43b) \d^L^v{t,x)\ < (Cic-t (C'feC'ie)2)(s/t)2-^'^s-3/2+'=^ 

When 0 ^ r/t ^ 3/5, we see that 4/5 ^ s/t ^ 1, then 

F{s) < (Cie+ (CfcCic)2)t-3/2sfe-5. 

Then, also by the sup-norm estimate (1X71) and (I7.15L we find that 

(7.44a) Idj^d^L-^vl < (Co.fe + + (CkCisfis/t)^-'^^ 3 -^^'^+'^^, 

(7.44b) Id^L'^vl < {Co,k + l)C'ie(s/t)^-'^‘^S-3/2 + fc5 

where we recall that Ci > Cq 

Then we conclude that there exists a positive constant C determined only by the structure of 
the system (11.111) such that 

(7.45a) 14.i(s) ^ C{Co,k + l)Cie3^^ + C{CkCiefs^\ 

(7.45b) 14.o(s) < C(Co,fe -b l)Ci£s'=^ + C{CkCiefs'^^. 

Now we consider together (17.421) and (17.451) and see that if Ck > 2C{Co,k + 1), then we can take 
ei := ■ Then we find that 

Vo.kis) ^ \CkCies>^^, Ih.fe(s) ^ Ukis) 

for all e ^ e/. This conclude that 32 ,k = S 2 so the case | J| = fc is proven. Then by induction 
we see that for all fc ^ N — A, (17.341) is established. Then taking £2 := min;s^A- 4 {efe} and 
Coo := niaxfc^Ar _4 Ck, we see that (I7.29ap and (I7.29bl) are established for all fc ^ — 4 and, more 

precisely, 

(7.46a) sup (t|L'^u|) ^ CooCi£s^^, 

Ms 

(7.46b) sup {{3/t)-^+'^^3^^‘^\dj_d^L-^v\) + sup ((s/t)-2+7-5s3/2|^/^J^|^ ^ CooCl£s'=^ 

Ms Ms 

(7.46c) sup {{s/t)-^+'^^3^^^\dad^L-’v\) ^ CooCi£s'=^ 

Ms 


From its definition, we see that Cqo is determined only from the structure of the system and 
therefore, we have proven (17.291) . □ 
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8. Refined energy estimate and completion of the bootstrap argument 


8.1. Overview. In this section, we derive the improved energy estimates (Ol) which concludes 
the main result. The improved estimates are classified in two categories. The first refers to the 
energy estimates of order higher than or equal to iV — 3, the second refers to those of order lower 
that or equal to — 4. 

First, we apply L'^ (with |/| + | J| ^ iV) to our system of equations 

(8.1) - ad^L'^u = L-’{do^vdpv) + Rd^L-^{v^), 

(8.2) - nd^L'^v + + H°‘^udcdpd^v = d^L'’]v. 

To be able to apply the energy estimate iProposition 12.11) . we need first to check Il2.12al) and 
(l2T^ . 

Lemma 8.1. There exists a positive constant Sq such that if the energy assumption (EH) is valid 
with Ci£ ^ 1 and e ^ £q, then the following estimates hold: 


(8.3) 

(8.4) 
with 


1 

— B < B < 2 B 

2 ^^ 771 ,c ^ ^ 9 ,c ^ 


r {s/t)\dcih'^^dtd^ L’^vdpd^ L‘’v\dx + f {s/t)\dth°‘^dad^ L'^vdpd^ L'^v\d: 


< M{s)E{s,d^L'^vY^‘^ 
M{s) sk 


j Cies-i/2+fc5^ - 3 ^ |J| + I J| ^ N, 

|ci£s-i+'=^ |/| + |J| ^ iV-4. 

Proof. The proof of (18.31) follows directly from (I7.17bl) . We remark that 

< Ci£t-^ < Ci£{s/tf 

where we have observed that ^ s ^ t in 3C. We get 

r \h°‘^dtd^L'’vdpd^L'’v\dx < Ci£ f \{s/tY‘dad^L'’vdpd^L'’v\dx < Ci£Eg^c{s,d^L'’v), 

r \h°'^dad^vdpd^L'’v\ < Ci£ f |(s/t)^(5o,(5^L'^P(5/3d'^L'^i;|da; < Ciei?g,c(s,d^L"^!;), 
J:Ks 

where we have used \{s/t)dad^L'^v\'^dx ^ Eg^ds, d^L'^v). 

So for some C" > 0 we have 

\Eg^c{s,d^L-^v) - Em,c{sd^L-^u)\ < C"C'i£i?g,c(s, d'^L'^u), 
and we choose Eq ^ 2 (yCi • Then, for e < eg, it holds 

|£'g,c(s,d^L'^u) - Em,c{s,S^L'’v)\ < CieEg^d^d^L'’v) iT;g,c(s,d^L'^u), 

which yields (18.3p . 

To derive (18.4p . we just need to observe that 


I 


|d.y/i“^dQ,d^L'^u|^da; 

(7i£ r f“^s“^(t/s)^ |(s/f)dad‘^L‘^u|^dx ^ Ci£ r ts“"^|(s/t)dQ-d‘^L'^u|^dx 


I'Ks 
„-2 ] 


< C'i£s ^£'g,c(s, d^L'^u), 
and we use the first estimate in (16.131) : 


L 




{C{Ci£fs-^+'^'^\ 7V-3^ |J| + |J| ^ iV, 

I C{Ci£) 


2„-2+2fc5 


|/| + |J| ^ lV-4. 
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So we see that 

r {s/t)\dah°‘^dtd^L'^vdpd^L-^v\dx ^ \\dah°‘'^dtd^L-^v\\L2^^)\\{s/t)di3d^ 

which is bounded by the right-hand side of (18.41) . The other term in the left-hand side is bounded 
in the same manner and we thus omit the details. □ 

8.2. Lower-order L? estimates. We remark that in lower order case where |/| 4-1 J| < iV —4, we 
have M{s) < and we need again the estimate on the source term L'^davdpv -I- 

Rv'^) and [H°‘^udadp,d^L'^]v. 

Lemma 8.2. Under the assumption of EH), the following estimates hold for |/| 4- | J| ^ iV — 4 
with I J| = k: 

(8.5) < (Cie)2s-3/2+'=^ 

(8.6) ||[Lr“^ud„d;3,d^L'']T||L2(Tc.) ^ {Ciefs-^+^^. 

Proof. The estimates of these terms relies on the basic and refined sup-norm estimates. We 
remark that 

^ “ Il+l2 = I ^ 

Jl + J2=J 

\l2\ + \J2\^N-5 

+ \\{t/s)dav\\^^^^J\{s/t)d^L'’di3v\\^2^^^-^ =: Ti 4- Ta- 

For Ta, we apply (|7.23bl) with |/| = 1 and (16. 8p and we conclude that 

Ta < Ci£(s/t)-i/2-^'^t-3/2Ci£sl'^l‘^£ < (C'i£)2s-3/2+fe5^ 

For Ti, we apply (|7.30p and (j6.9p and we conclude that 

Ti < C'i£(s/t)"l/2"^'^t"3/2sl'^il'^Ci£sl'^"l^£ < (C'i£)^s"3/2 + '=^ 

The estimate on the term L'^ (v'^) is similar by apply (I7.29bl) (16.811 and we omit the details. 
To see the estimate on [H^^^udadp, L'^]v is quite similar, we just need to remark that it is a 
linear combination of the following terms: 

L'^^ud^ L'^^dadpv, d^^L'^^ud^^L'^^dadpv, udadpd^L'^^v 

where Ii + h = I, Ji 4- Ja = T, 4- J 2 = J with \ J[ \ > 1, |/i| > 1 and | J 2 | < | J| — 1. For the last 
term we apply (17.17bl) and (16.31) : 

Indeedik\\{t/s)u\\^^^^^^^\\{s/t)do,dpd^L-’'iv\\^2^^^^ 

<Cies-^Cies^^ < {Cief . 

For the first term, we see that | J(| ^ TV — 4, then we apply (I7.29ap and (16.3p : 

<C'i£s-i+l'^il^C'i£sl‘^=l'^ < {Cief 3-^+’^^. 

For the second term, we see that when |/i| = 1 and Ji = 0, 

Wd^ud^'^ L-’ dadi3v\\L2(^jl^) ^||(t/s)d^M||Lco(M^) \\{s/t)d^ L-’ dadpv\\L 2 (j^^-.^ 

<Ci£||(t/s)t“^/^S“^||icx,(5^^) C'i£s'''^ ^ (C'i£)^s“^ + '=^. 
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When |/i| + I Ji| ^ 2, we see that I/ 2 I + | J 2 I ^ iV — 6. Then by the first inequality in (16.lip and 
(lOD . we find 

and we conclude with (18.61) . □ 

8.3. Higher-order estimates. When N — 3 ^ |/| + | J| ^ N, the energy estimate is more 
complicated. For the source terms we have the following estimates. 

Lemma 8.3. Under the energy assumption dSU) the following estimates hold for N—4 ^ |/| + | J| ^ 
N and | J| = k: 

(8.7) \\d^LJ{P^f^dc,vdpv)\\^,^^^^ + ^ , 

(8.8) < (Cie)2s-i/2+fc-5. 

Proof. The proof relies on the refined decay estimate (17.291) and the basic estimates. We begin 
with (18.71) . We remark that L'^[davdjsv) is a linear combination of the following terms 

with h + I 2 = I, Ji + J 2 = J- We see that when |/i| + | Ji| = 0, we apply (I7.23bl) on daV ( with 
1 ^ TV — 4) and (16.51) 

\\d^^L‘^^davd^^L‘^^d/ 3 v\\^ 2 ^^^^ = \\{t/s)dav{s/t)d^ P-^ 

< C'ie||(t/s)(s/t)^/^“'^‘^t“^/^||ioo(3<:«)C'C'ies^/^+''^ 

< {Ciefs-^+'^^. 

When 1 ^ |Ji| + |Ji| ^ iV —4, we see that 4 ^ I/ 2 I + IJ 2 I ^ iV — 1. Then we apply (I7.29cl) and the 
third inequality in (16.6F 

<C'fcC'ie||(s/t)-l/2-75^-3/2+|Ji|5||^_^^^^^(^^^gl/2+|J2|5 

<Ck{Ciefs-^+>^^. 

When iV — 3 ^ |/i| + I Ji| < — 1, we see that 1 ^ |72| + | J 2 I ^ 3 ^ — 4. Then we apply the 

third inequality in (16.61) and (l7.29cF Similar to the former case, 

When |/i| + |Ji| = N and I/ 2 I + IJ 2 I = 0, the estimate is derived similarly as in the first case by 
exchanging the role of daV and djsv. The we conclude that 

< (C7ie)2s-i+'='5. 

The estimate on d^L'^(u^) is quite similar by applying (16.51) and (I7.29bl) . we omit the detail. 

The estimate on [H°^l^udadp,d^L'']v is as follows: we observe that this term is a linear combi¬ 
nation of the following terms 

L''^ud^L'’'^dadpv, d^^L'’^ud^'^L'’‘^dadpv, udadpd^^L'^'^ v 

where Ii + h = I, Ji + J 2 = J, J'l + J 2 = J with \ J[\ ^ 1 |/i| 5= 1 and | J 2 | < | J| — 1. The last 
term is bounded by applying (17.171) and (16.51) : 

<C'i£s-iC'iesi/2+l'^=l^ < (C7ie)^s-^/2+'=‘^. 
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For the first terra, we make the following observation. When 1 ^ \J[\ ^ TV — 4, we have 
4 ^ |/| + IJ 2 I ^ fV — 1. Then we apply (j7.29al) and (16.71) : 

<C'i£s"^sl'^f^Cies^/2+l“'2l'^ < {Cief 3 -^!'^+^^. 

When N — 3 ^ | J{ | ^ N, we see that |/| + |J 2 | ^3^iV — 5. Then we apply the Hardy inequality 
in the form (I6.20a|) as well as (17.311) . So we see that 

The second term is easier, since the factor has better decay when |/i| 5= 1. Then we see 

that when |/i| = 1 and | J 2 I = 0, 

Wd^^ud^'^L'^dadpv\\L-2(^^^) ^ \\{t/s)d^^u{s/t)d^'^L-^ 

when 2 ^ |/i| + | Ji| ^ N — 2, I/ 2 I + |•/ 2 | ^ N — 2. Then we apply the third inequality in (16.6p and 
we see that 

When TV — 1 ^ |/i| + I Ji| ^ TV, I/ 2 I + | J 2 I sS 1 < TV — 7 then we apply (16.51) and (I7.33F Then, we 
obtain 

which completes the argument. □ 


8.4. Proof of Proposition [5711 Our aim is to establish the improved energy estimate (j5.2p and 
to conclude the proof of Theorem ll.il that is, we now establish Proposition 15.11 The strategy is 
to apply the energy estimate 12.II with (18.31) . (18.4p . (18.51) . (18.61) . (18.71) . and (18.8F 

We need to specify the constants and we denote by C a sufficiently large constant determined 
only by the structure of the system such that all of the above estimates hold true. We derive the 
wave equation of (jl.lip hy d^L'^: 

-nd^L-^u = d^L-’{P°‘^dc.vdi3v) + d^L-^{v^). 

Recall the energy estimate (12.101) 


Em{s,d^L-’uYl^ ^ Em{2,d^L^u 


'■""f 


ll□w||L2(M^) ds 


with ||□M||L 2 (J^_) ^ ||(5-'^L'^(P“^dQ,'yd;3u)||j^2(j^^) + \\d ^Then by (|531) . when |/| + 
I J| ^ TV — 4, we have ||□^^||L 2 (M^) ^ C(C'i£)^s“^/^+^^, and we conclude that 


(8.9) 


P„^(s, d^L^uf!^ ^ CC'o£ + C'(C'i£)^ 


When TV — 3 ^ |/| + | J| ^ TV and | J| = fc, by (18.71) 


( 8 . 10 ) 


P™(s, d^L^uf/'^ ^ CCos + C{Cief J 5 "^+'='^ ds 
^CCoe + C{Ciefs'‘^. 


For the energy estimates on v, we apply to the Klein-Gordon equation in (11.111) and obtain 
—n\d^L'^v + ]^Jy _|_ c^d^L'^v = [E[°'dudadp,d^ L^^v. 
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Then by (|2.13|) . and (18.31) (with k = 1/2), we find 

Em,c{s, L '^^ L '^^ E ^ I|-^(®)IIl^(Ms) ds. 

When |/| + |J|^7V — 4, we rely (18.61) and (18.41) and observe that 

E^^a{s,d^L^vy/‘^ ^CCos+C{Cief ^ 5 -^+'=^ ds 
^CCoe+C{Ciefs'^^. 

When TV — 3 ^ |/| + I J| ^ N, we apply (18.81) and (18.41) and observe that 
E,n c(s, ^ CCos + C{Cief f 5 - 1 / 2 +^^ 

( 8 . 12 ) J2 

iiCCoe + C{Ciefs^^^+'^^. 

Finally, by choosing Ci > ACCq and e ^ (dCCi)"^, (18.91) -- (18. 121) lead to (15.2F 
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